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Abstract 

The paper presents the geometry of Lie algebroids and its apph- 
cations to optimal control. The first part deals with the theory of 
Lie algebroids, connections on Lie algebroids and dynamical systems 
defined on Lie algebroids (mainly Lagrangian and Hamiltonian sys- 
tems). In the second part we use the framework of Lie algebroids in 
the study of distributional systems (drift less control affine systems) 
with holonomic or nonholonomic distributions. 
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PREFACE 



The present paper is devoted to Lie algebroids geometry and its ap- 
plications to optimal control and variational calculus. The framework of 
the differential geometry is very useful in modelling and understanding of a 
large class of natural phenomena. The Lie geometric methods are applied 
successfully in differential equations, optimal control theory or theoretical 
physics. In the most of cases the study is starting with a variational problem 
formulated for a regular Lagrangian (see [1]), on the tangent bundle TM 
over the manifold M and very often the whole set of problems is transferred 
on the dual space T*M, endowed with a Hamiltonian function, via Legen- 
dre transformation. The case of a non-regular Lagrangians is also studied. 
The problem in this case is that the proposed Lagrangian formalism yields 
a singular Lagrangian description, which makes the Legendre transform ill- 
defined and thus no straightforward Hamiltonian formulation can be related. 
One of the motivations for the present work is the study of Lagrangian sys- 
tems subjected to external constraints (holonomic or nonholonomic) . These 
systems have a wide application in many different areas as optimal control 
theory, mathematical economics or sub-Riemannian geometry. 

In the last years the investigations have led to a geometric framework 
which is covering these phenomena. It is precisely the underlying structure 
of a Lie algebroid on the phase space which allows a unified treatment. 
This idea was first introduced by A. Weinstein [1251 [23] in order to define 
a Lagrangian formalism which is very useful for the various types of such 
systems. 

The concept of Lie algebroids have been introduced into differential ge- 
ometry since the early 1950, and also can be found in physics and algebra, 
under a wide variety of names. However, the fundamental concept has been 
introduced in sixties by J. Pradines |112j in relation with Lie groupoids. For 
every Lie groupoid there exists an associated Lie algebroid, like as for every 
Lie group there exists an associated Lie algebra. A Lie algebroid [73l [75] 
over a smooth manifold M is a real vector bundle {E, vr, M) with a Lie alge- 
bra structure on its space of sections, and an application a, named anchor, 
which induces a Lie algebra homomorphism from sections of E to vector 
fields on M. It is convenient to think a Lie algebroid as a substituent for 
the tangent bundle of M, an element e of as a generalized velocity, and 
the actual velocity v on TM is obtained when applying the anchor to e, i.e., 
a{e) = V. 

The basic example of Lie algebroid over the manifold M is the tangent 
bundle TM itself, with the identity mapping as anchor. Every integrable 
distribution of TM is a Lie algebroid with the inclusion as anchor and in- 
duced Lie bracket, and every Lie algebra is a Lie algebroid over one point. 
An important Lie algebroid is the cotangent bundle of a Poisson manifold 
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|60j . Being related to many areas of geometry, as connections theory [75\ HOl 
[Ml [3 H H9 l [76 l [TTH [95] cohomology [73l[75] foliations and pseudogroups, sym- 
plectic and Poisson geometry [Ml UMl [Ml [Ml [HD [SH [Ml IMl IM IM [W] 
the Lie algebroids are today the object of extensive studies. More precisely, 
Lie algebroids have applications in mechanical systems and optimal control 
theory [Ml[ZZl[ZHl[l9l[29l[lll[l[92l[96l[98l[l6] (distributional systems) and 
are a natural framework in which one can be developed the theory of dif- 
ferential operators (exterior derivative and Lie derivative) and differential 
equations. 

In his papers [73l [75] K. Mackenzie has been achieved a unitary study 
of Lie groupoids and algebroids and together with P. Higgins [36] have in- 
troduced the notion of prolongation of a Lie algebroid over a smooth map, 
useful in the study of induced vector bundle by the Lie algebroid structure. 
Using the geometry of Lie algebroids, A. Weinstein |125] shows that is possi- 
ble to give a common description of the most interesting classical mechanical 
systems. He developed a generalized theory of Lagrangian mechanics and 
obtained the equations of motions, using the Poisson structure on the dual 
of a Lie algebroid and Legendre transformation associated with a regular 
Lagrangian. In the last years the problems raised by A. Weinstein have 
been investigated by many authors. Thus, E. Martinez [BS1[7DI[72] obtained 
the same Euler-Lagrange equations using the symplectic formalism for La- 
grangian and Hamiltonian, similarly with the J. Klein formalism |57] for the 
classical Lagrangian mechanics. 

In the classical version of the tangent bundle {E = TM) the Klein's 
method is based on the vector bundle structure of TAf and the existence 
of a vector-valued 1-form. Such a form does not exist for a general Lie al- 
gebroid [65] because of different dimensions of the horizontal and vertical 
distributions, and so Klein's approach is not applied directly. To overcome 
this difficulty, E. Martinez, M de Leon, J. C. Marero [MllSl] have proposed a 
modified version, in which the bundles tangent to E and E* are replaced by 
the prolongations TE and TE* (in sense Higgins and Mackenzie @S])- The 
nonholonomic Lagrangian systems and Hamiltonian mechanics on Lie alge- 
broids are studied by a group of E. Martinez [64j. The first step in studying 
the mechanical control systems on Lie algebroids seems to be done by J. 
Cortes and E. Martinez [30], which also approached the problem of accessi- 
bility and controllability. A framework for nonholonomic systems, using a 
subbundle of a Lie algebroids is proposed by T. Mestdag and B. Langerock 
[77] . A start in the study of some problems of control affine systems and 
sub-Riemannian geometry, using the framework of Lie algebroids is due to 
D. Hrimiuc and L. Popescu [Ml [Ml [SB]- 

Control theory is splitting in two major branches: the first is the control 
theory of problems described by partial differential equations where the ob- 
jective functionals are mostly quadratic forms, and the second is the control 
theory of problems described by the parameter dependent ordinary differ- 
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ential equations. In this last case it is more frequent to deal with non-linear 
systems and non-quadratic objective functional. The mathematical models 
from the optimal control theory cover also the economic growth in both open 
and closed economies, exploitation of (non-) renewable resources, pollution 
control, behavior of firms or differential games \38 \ 11151 1116j . 

The geometric methods in the control theory have been applied by many 
authors (see [THl EUJ [TSl EI])- One of the most important issues in the 
geometric approach is the analysis of the solution to the optimal control 
problem as provided by Pontryagin's Maximum Principle; that is, the curve 
c{t) = {x{t),u{t)) is an optimal trajectory if there exists a lifting of x{t) to 
the dual space {x{t),p{t)) satisfying the Hamilton equations, together with 
a maximization condition for the Hamiltonian with respect to the control 
variables u{t). 

In the paper [7l] E. Martinez presents the Pontryagin Maximum Princi- 
ple on Lie algebroids using the prolongation (in sense of Higgins and Macken- 
zie [l6]) of the Lie algebroid over the vector bundle projection of a dual 
bundle. In this paper we study some distributional systems with positive 
homogeneous cost, using the Pontryagin Maximum Principle at the level of 
a Lie algebroid. 

"In spite of that, the control theory can be considered part of the gen- 
eral theory of differential equations, the problems that inspires it and some 
of the results obtained so far, have configured a theory with a strong and 
definite personality, that is already offering interesting returns to its an- 
cestors. For instance, the geometrization of non- linear affine-input con- 
trol theory problems by introducing Lie-geometrical methods into its anal- 
ysis, started already by R, Brocket [18], is inspiring classical Riemannian 
geometry and creating what is called today sub-Riemannian geometry" 

[nHiiHaiiiEiiiiiiiiiiniiiiiiiaiis]. 

If M is a smooth n-dimensional manifold then a sub-Riemannian struc- 
ture on M is a pair {D,g) where D is a distribution of rank m and (7 is a 
Riemannian metric on D. A piecewise smooth curve on M is called hori- 
zontal if its tangent vectors are in D. The length of a horizontal curve c is 
defined by 

L{c) = ^Mm)dt, (1) 

where g is a Riemannian metric on D. The distance between two points a 
and b is d{a,b) = infL{c), where the infimum is taken over all horizontal 
curves connecting a to b. The distance is assumed to be infinite if there is no 
horizontal curve that connects these two points. If locally, the distribution 
D of rank m is generated by Xj, z = 1, m a sub-Riemannian structure on M 
is locally given by a control system 

m 

x = J2Mt)Xi{x), (2) 
1=1 
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of constant rank m, with the controls u{.). The controlled paths are obtained 
by integrating the system (2) and are the geodesies in the framework of 
sub-Riemannian geometry. If D is assumed to be bracket generating, i.e. 
sections of D and iterated brackets span the entire tangent space TM, by a 
well-known theorem of Chow |25j the system (2) is controllable, that is for 
any two points a and b, there exists a horizontal curve which connects these 
points (M is assumed to be connected). 

The concept of sub-Riemannian geometry can be extended to a more 
general setting, |49t f27\ [2H] by replacing the Riemannian metric with a posi- 
tive homogeneous one. For the theory of optimal control this extension 
is equivalent to the change of the quadratic cost of a control affine system 
with a positive homogeneous cost. Also, the case of distribution D with non- 
constant rank is generating interesting examples (Grushin case \37\ I49j). 

The case when the distribution D generated by vector fields Xi, i = 
l,m is integrable is also studied. In this case the distribution determines a 
foliation on M and two points can be joined if and only if belongs to the 
same leaf. In order to find the optimal trajectory of the system one uses the 
Pontryagin Maximum Principle at the level of Lie algebroids, built different 
in the case of holonomic or nonholonomic distributions. 

* * * 
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The paper is organized in two parts. The first part entitled The ge- 
ometry of Lie algebroids contains seven chapters. In the first chapter some 
prehminaries concerning geometrical structures on the total space of a vector 
bundle are presented [82]. We focus on the notions of nonlinear connection 
and covariant derivative. In the next chapter we present the notion of Lie 
algehroid including the cohomology and structure equations (73]. The notion 
of prolongation of a Lie algebroid over the vector bundle projection is stud- 
ied in the chapter three. The Ehresmann nonlinear connection N = —CsJ 
with the coefficients given by 



is investigated and the relations with the Ehresmann connections on tangent 
bundles TE and TM are pointed out. In the chapter four we introduce the 
notion of dynamical covariant derivative and metric nonlinear connection 
at the level of the Lie algebroid TE. The Lagrangian formalism on Lie 
algebroids yields a canonical semispray [70] 



and a canonical Ehresman connection, which is a metric nonlinear connec- 
tion. We also have the Lagrange equations on Lie algebroids given by |125j 



In the case of positive homogeneous Lagrangian (Finsler function) we find 
a canonical Ehresmann connection which depends only on Finsler function 
and the structure functions of the Lie algebroid. 

In the chapter five we deal with the prolongation of a Lie algebroid over 
the vector bundle projections of a dual bundle. We introduce the notions of 
dual adapted tangent structure J and J'-regular sections. These structures 
induce a canonical nonlinear connection M = —LpJ with the coefficients 
given by [50] 



In the case of Hamiltonian formalism these coefficients become [103] 
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where {•,•} is the Poisson bracket. The corresponding Hamilton equations 
on Lie algebroid are given by \125\ [M] 



'dt 



j dV. dfi, 



,i&H_ 




dn 



a 



dfj,, 



'a 



dt 



In the chapter six we introduce the notion of dynamical covariant derivative 
and metric nonlinear connection at the level of a Lie algebroid TE*. We 
prove that the canonical nonlinear connection induces by a regular Hamilto- 
nian is a unique metric and symmetric nonlinear connection. In the chapter 
seven we investigate some aspects of the Lie algebroids geometry endowed 
with a Poisson structures, the so-called Poisson-Lie algebroids. 

Author's papers [ia[50l[9ll[92l[9ll95l[96l[97l[99l [TOOl [Ml [TOll [T05l 
llU6t \107\ llU9j are used in writting this part. 

The purpose of the second part entitled Optimal Control is to study 
the drift less control affine systems (distributional systems) with positive 
homogeneous cost, using the Pontryagin Maximum Principle at the level of 
a Lie algebroid in the case of constant rank of distribution. 

We prove that the framework of Lie algebroids is better than cotangent 
bundle in order to solve some problems of drift less control affine systems. 
In the first chapter the known results on the optimal control systems are 
recalled by geometric viewpoint. In the next chapter the distributional sys- 
tems are presented and the relation between the Hamiltonians on E* and 
T*M is given by 



We investigate the cases of holonomic and nonholonomic distributions with 
constant rank. In the holonomic case, we will consider the Lie algebroid 
being just the distribution whereas in the nonholonomic case (i.e., strong 
bracket generating distribution) the Lie algebroid is the tangent bundle with 
the basis given by vectors of distribution completed by the first Lie brackets. 
Also, the case of distribution D with non-constant rank is studied in the 
last two sections of the chapter and some interesting examples are given. In 
the last chapter we present the intrinsic relation between the distributional 
systems and sub-Riemannian geometry. Thus, the optimal trajectory of our 
distributional systems are the geodesies in the framework of sub-Riemannian 
geometry. We investigate two classical cases: Grusin plan and Heisenberg 
group, but equipped with positive homogeneous costs (Randers metric). We 
are using the Pontryagin Maximum Principle at the level of Lie algebroids, 
in the case of Heisenberg group and show that this idea is very useful in 
order to solve a large class of distributional systems. Author's papers [501 
[Ml [Ml Unil [Ml IMl IM] are used in writing this part. 

In my opinion, the paper is useful to a large class of readers: graduate 
students, mathematicians and to everybody else interested in the subject 



H{p)='H{fi), f^ = a^p), peT^M, fi^E*. 



9 



of differential geometry, differential equations or optimal control. I want to 
address my thanks to all authors mentioned in this paper and to everybody 
else I forgot to mention, without any intention, in the Bibliography. 

Finally, I wish to address my thanks to the referees for many useful 
remarks and suggestions concerning this paper. I should like to express 
the deep gratitude to professor D. Hrimiuc for the collaboration during the 
postdoctoral fellowship at the University of Alberta, Edmonton, Canada, 
where many ideas presented in this paper have been started. Also, I want 
to address my thanks to Professor P. Stavre for support and guidance given 
me in life and in mathematics. 

Acknowledgments: This work was supported by the strategic grant POS- 
DRU/89/1.5/S/61968, Project ID61968 (2009), co-financed by the European 
Social Fund within the Sectorial Operational Program Human Resources 
Development 2007-2013. 

Craiova, December 2012 

Assoc. Prof. Liviu Popescu 
University of Craiova 
Department of Applied Mathematics 
13 "Al. I. Cuza" St., Craiova 200585 
e-mail: liviupopescu@central.ucv.ro 



10 



1 THE GEOMETRY OF LIE ALGEBROIDS 



The purpose of this first part is to study the geometry of a Lie algebroid 
and its prolongations over the vector bundles projections. A Lie algebroid 
[73\ [75] over a smooth manifold M is a real vector bundle {E, it, M) with a 
Lie algebra structure on its space of sections, and an application cr, named 
the anchor, which induces a Lie algebra homomorphism from the sections 
of E to vector fields on M. For this reason, in the first chapter we present 
some results on the geometry of the total space of a vector bundle, including 
nonlinear connections and covariant derivatives. In the next chapter we give 
only the relevant formulas for Lie algebroid cohomology we shall need later, 
and refer the reader to the monograph [73j for further details. 

The chapter three deals with the prolongation TE of a Lie algebroid 
over the vector bundle projection. We introduce the Ehresmann nonlinear 
connection on the Lie algebroid TE and study its properties [95l [92] . We 
show that the vertical part of the Lie brackets of horizontal sections from 
the basis represents the components of the curvature tensor of the nonlinear 
connection. We study the related connections and show that a connection 
on the tangent bundle TE induces a connection on the Lie algebroid TE. 
We introduce an almost complex structure on Lie algebroids and prove that 
its integrability is characterized by zero torsion and curvature property of 
the connection. We present the notion of dynamical covariant derivative 
at the level of a Lie algebroid and show that the metric compatibility of 
the semispray and associated nonlinear connection gives the one of the so 
called Helmholtz conditions of the inverse problem of Lagrangian Mechanics. 
In the homogeneous case a canonical nonlinear connection associated to a 
Finsler function is determined. We study the linear connections on TE and 
determine the torsion and curvature. 

In the chapter four we study the dynamical covariant derivative and 
metric nonlinear connection on TE [104] . We introduce the dynamical co- 
variant derivative as a tensor derivation and study the compatibility condi- 
tions with a pseudo-Riemannian metric. In the case of SODE connection we 
find the expression of Jacobi endomorphism and its relation with curvature 
tensor. We prove that the canonical nonlinear connection induced by a reg- 
ular Lagrangian is a unique connection which is metric and compatible with 
symplectic structure. Also the invariant form of the Helmholtz conditions 
on Lie algebroids are given (see also [35]). 

The chapter five deals with the prolongation TE* of a Lie algebroid over 
the vector bundle projection of a dual bundle. We study the properties of the 
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connections on TE* [50\ [Ml [96] and introduce the notions of adapted almost 
tangent structure, almost complex structure and characterize the integrabil- 
ity conditions in terms of torsion and curvature of the connection. We prove 
that every ^-regular section (in particular, any regular Hamiltonian on E*) 
determines a canonical Ehresmann connection on the Lie algebroid TE*. 
We introduce some generalizations of the Hamilton sections, as a mechani- 
cal structures and semi-Hamiltonian sections and study their properties. In 
the last part of this chapter, using the diffeomorphism from TE* and TE 
|64| induced by a regular Hamiltonian, we can transfer many geometrical 
results between these spaces. Thus, a semispray on TE is transformed into 
a semi-Hamiltonian section on TE* if and only if the nonlinear connection 
on TE determined by semispray is just the canonical nonlinear connection 
induced by regular Lagrangian, via Legendre transformation. 

In the chapter six we study the dynamical covariant derivative and met- 
ric nonlinear connection on TE* [109j . Using the notion of ^-regular section 
we introduce the dynamical covariant derivative as a tensor derivation. In 
the case of nonlinear connection induce by a j7-regular section we find the 
expression of Jacobi endomorphism and its relation with curvature tensor. 
Finally, we prove that the canonical nonlinear connection induced by a reg- 
ular Hamiltonian is the unique metric and symmetric nonlinear connection. 

In the chapter seven we investigate some aspects of the Lie algebroids 
geometry endowed with a Poisson structures [971 EHl I107j , which generalize 
the Poisson manifolds. We recall the Cartan calculus and the Schouten- 
Nijenhuis bracket at the level of Lie algebroids and introduce the Poisson 
structure on Lie algebroids. We study the properties of linear contravariant 
connection and its tensors of torsion and curvature. In the last part of this 
section we find a Poisson connection which depends only on the Poisson 
bivector and structural functions of Lie algebroid, which generalize some 
results of Fernandes from [3S]. Also the geodesic equations are given. We 
study the properties of the complete lift of a Poisson bivector on TE and 
introduce the notion of horizontal lift. The compatibility conditions of these 
bivectors are investigated. Finally, the compatibility conditions between the 
canonical Poisson structure and the horizontal lift on TE* are given. 
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1.1 The geometry of the total space of a vector bundle 
1.1.1 Connections on vector bundles 

The connections theory is an important topic of the differential geometry 
with important applications in Differential Equations or Optimal Control. 
In this section we shall present only the notion of Ehresmann nonlinear con- 
nection and induced geometrical structures. Also, the covariant derivative is 
described. For more details and complete proofs we refer to the monographs 
[591 Ea [791182]. 

Let us consider the n-dimensional differentiable manifold M and a vector 
bundle {E,Tr,M) over M, with type fibre F = W^. The structure of the 
vector bundle is given by a vectorial atlas {(C/j, ^j, such that 

1) {Ui)i^i is an open covering of the manifold M. 

2) The mappings -i/'j : '^~'^{Ui) Ui x are bijective and satisfy the 
relation 

7r{i^-\xJ)) =x, :eGM, 

3) For every pair £ I x I, such that Ui n C/j 7^ there exists a 
smooth mapping gij : Ui D Uj — GL(m,M) with = ^'J^ o gji{x) for 
every x € UiD Uj where is the restriction of to {x} x R"^. 

Let {{Ui,ipi)}ii^j be an atlas on the manifold M such that Ui belongs 
to the maps domain into vectorial atlas {{Ui,ipi,M.'^)}iQj . We obtain that 
{(7r"^(C/i),/ii)}ig/ with 

h, : 7r-\Ui) ^ M" X M™, hi{u) = (v9,(7r(^.)), (^i)), 

is a differentiable atlas on the manifold E. If (Uj, ipj) is another local chart 
on M with U nUj ^9 and {Uj,'ipj,R"') is a bundle chart, then 

{hj o h-^){x,y) = {{ifj o ipi){x),gji{ip^{x)y) , (x,y) G x M*". 

We will denote by (x) the entries of the matrix associated with the linear 
application gji{ipi{x)). For x £ U C M we take ifi{x) = (x*) G M", i = l,n 
and considering (pj{x) = (x* ) then (fj o ipT^ has the form 

= X* (x^, x"), rank i = n. 
^ ' ' ^' ydx' j 

A tangent vector at the point x € M will be locally represented by the 
pair (x*,?/*), where (y*) is given by Xx = U^-^- Therefore, the transforma- 
tions of coordinates on differentiable manifold TM have the form 

X* = X* (x\...,x"), "rankl-^] =n, = -Q-fV'- (l-l-l) 
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A point u £ E on the total space of a vector bundle {E, vr, M) is locally rep- 
resented by the pair (x*, y") G M"' x . The transformations of coordinates 
{x^, y"") — )• (x* , ) on are of the form [82] 

X*' = x*'(x^, x"), ranA: ( ) = JT^, 

(1.1.2) 

y'*' = M^'(x)y'^, rank{M^' (x)) = m. 

For the vector bundle (i?, vr, M) we consider vr* : TE — t- TM the tangent 
application of vr. The application vr,,, is a vr -morphism of vector bun- 
dles between tangent bundles {TE^tte, E) and {TM,-km, M). The kernel 
of this TT-morphism is a subbundle of vector bundle {TE^tte, E), denoted 
(yE,iTv,E) = Ker vr*, which will be called the vertical subbundle. The 
total space is VE = UVu , where Vu = Ker tt^, u £ E. 

A tangent vector Xu at the point u £ E has the local representation 
{x\ y'^,X\A''), where the coefficients (X*) G and {A") E R™ are defined 
by the relation Xu = + ^"afs"- The tangent application vr* is locally 

represented by vr*(x,y,X, ^) = {x,X). Consequently, the local fibres of 
vector bundle {TE, vr* , TM) are isomorph with {x} x R™ x {X} x R™ ~ R^'". 
Because vr*(g|s-) = 0, it results that the functions {-^}, a = l,m determine 
a local basis of the vertical distribution {u ^ Vu \ u £ E}, which means 
that Vu is integrable. The elements of the vertical subbundle have de form 
{x,y,0, A), that is the fibres of vertical subbundle VE are locally isomorph 
with R™. Let ■n*TM be the induced vector bundle of tangent bundle over 
the application tt : E ^ M and vr! : TE vr*rM given by 

vr!(X„) = (n,vr*(X„)). 

This application is a morphism of vector bundles, and follows that the ap- 
plication vr! is a surjection and 

Ker vr! = Ker vr* = VE. 

Therefore, it results that the following sequence of vector bundles over E is 
exact 

O^VE Ate ^7T*TM ^0. (1.1.3) 

where i : VE — )• TE is the inclusion map. 

We can present now a definition of the Ehresmann connection, called 
usually nonlinear connection. 

Definition 1.1.1 The nonlinear connection in the vector bundle {E,tt,M) 
is a splitting on the left of the exact sequences (1.1.3). 

It results that a nonlinear connection in {E, vr, M) is a morphism of 
vector bundle N : TE — )• VE such that N oi = Id\vE- The kernel of the 
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morphism is a vector subbundle of the bundle {TE, tte, E), which will be 
called horizontal subbundle and will be denoted {HE, tth^E). Therefore, the 
vector bundle {TE, tte^E) is the Whitney sum of the horizontal and vertical 
subbundles. Thus we have the following characterization of the nonlinear 
connection. 

Proposition 1.1.1 A nonlinear connection in {E,Tr,M) is determined by 
the existence of the vector subbundle {HE,tth,E) of the tangent bundle over 
E, (TE, TTE,E) such that TE = HE VE. 

The restriction vr! {he of the application vr! to the horizontal subbundle 
HE is an isomorphism of vector bundles. The component vr* : HE — t- TM 
of the application vr! \he is a tt- morphism and its restriction to the fibres 
is an isomorphism. Therefore, for any vector field X on M, there exists a 
horizontal vector field on E, such that tt^:{X^) = X. The vector field X^ is 
called the horizontal lift of the vector field X. The horizontal lift has a local 
representation 

d 6 

I = l,n. 



thus, we determine a local basis of HuE. This vector field can be 

represented in the form ^ = + Bf^, but the condition vr* (^) = 

^ implies A^- = 5*- and the fact that ^ are the kernel of the mapping 
gives Bf = —Nf-. It results 

A = A_Arf^A (114) 
5xi dx^ ' dy^' ^ ' 

We obtain a new basis gf^) tangent bundle T^E which is called 

the Berwald basis associated to the nonlinear connection N. 

If we denote N^, i = l,n, a = l,m the coefficients of a nonlinear con- 



nection, then it results 

Proposition 1.1.2 By a change of the local coordinates (1.1.2) on the vec- 
tor bundle {E,tt,M), the local coefficients N?- of a nonlinear connection N 
change as follows 



^^^ = ^«"^^-^^'^- 

A nonlinear connection in the vector bundle (£'\{0}, vr, M) is said to be 
homogeneous (respective linear) if the coefficients N^{x, y) are homogeneous 
(respective linear) with respect to the second argument. 

We define a morphism v : X{E) — > X{E) such that v{X) = —N{X) if 
X e T{VE) and v{X) = 0, for X e T{HE). It follows that the morphism 
V : X{E) — )• X{E) has the properties: 
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1° v(X{E)) c r{VE), 

2° {v{X) = X}^X £ T{VE), 
and it determines a nonlinear connection in (E, ir, M) since Kerv = T{HE). 

It results that v"^ = v, i.e. u is a projector, which is called the vertical 
projector of the nonlinear connection N. Analogously, a nonlinear con- 
nection in the vector bundle {E, vr, M) is characterized by the morphism 
h : X{E) X{E) with the properties h'^ = h, Ker h = TiVE), and h is 
called the horizontal projector of the nonlinear connection A^. It results that 
h + V = Id. In [82] one proves: 

Proposition 1.1.3 A nonlinear connection in a vector bundle {E,tt,M) 
is characterized by an almost product structure P on E whose distribution 
of eigensub spaces which correspond to the eigenvalue —1 coincides to the 
vertical distribution. 

Prom the previous consideration it follows that 

P = 2h- Id = Id-2v = h-v,P'^ = Id. 

The existence of the nonlinear connection in [E, vr, M) leads to the decom- 
position 

X{E) = T{HE) e T{VE) => X = hX + vX, VX € X{E). 
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1.1.2 Covariant derivative 



We consider a local basis {sa} of the sections of the bundle 7r^^(f7) — t- U, 
U C M and it results p = p"'{x)sa, p & r(£'). The nonlinear connection 
induces a covariant derivative for the sections in E defined as follows 

Dxp = X^ (^^ + Nt{x,pix))X' 
where 

X = X'^^eTiTM). 

This covariant derivative has the properties: 

1° Dx+Y = Dx + Dy, 

2°Dfxp = fDxP, X,Yer{TM). 
If N is homogeneous one also has 

3° Dxifp) = X{f)p + fDxP, p G nE), f e F{M). 
If N is linear, then the covariant derivative satisfies 1°, 2°, 3° and 

4° Dx{pi+P2) = Dxpi + Dxp2, Pi,P2 e r(^). 

Let c : [a, b] — )• M, t c(t) be a smooth curve on the manifold M and 
c: [a, b] — >■ TM the tangent field along the curve c. Setting D-^p = we say 
that the section p in {E, vr, M) is parallel along c if ^ = and it results: 

Proposition 1.1.4 A local section p of the bundle 7r~^{U) U , U d M 
is parallel along the curve c with c{[a,b]) C U if and only if 

do"' dx^ 

-^+N^{x,pix))—=0, p = p''{x)sa, xeU. 

The decomposition T^E = H^E © VuE in the vector bundle {E, tt, M) 
leads to the decomposition {T^E)* = {V^E)^ © {HuE)^ , for u G E, where 
{VuE)-^ denotes the subspace of 1-forms on T^E which vanish on horizontal 
vectors, and [HuE)-^ denotes the subspace of 1-forms on T^E which vanish 
on vertical vectors. Therefore, we have 



huj + vuj, yu£A^{E), 



where hu){X) = oj{hX) and vijj{X) = ijj{vX). The dual Berwald basis is 
(dx*, where 

5y« = dy" + N^dx\ (1.1.6) 

Proposition 1.1.5 The Berwald basis (^-^i-^^ o,nd its dual {dx^,6y°') 
transform under a change of coordinates (1.1.2) as follows 

_S_^dx^_5_ d ^ , d 
5x^ dx' 6x'' ' " dy""' ' 
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dx^ 
dx^' 



dx' 



5y- = M^^Sy'^ . 



The Lie brackets of the vector fields of Berwald basis are given by 



_6_ _6_ 

6x^ ' 6x^ 



_6_ _d_ 

6x^ ' dy'^ 



where 



dN^ d 

Qya Qyb ' 

5m 5Nt 



d _d_ 

Qya ' Qyh 



6x^ 



= 0, 
(1.1.7) 

(1.1.8) 



In the following we deal with the curvature of a nonlinear connection N. 

Definition 1.1.2 The curvature of a nonlinear connection N is given by 

= -N„, (1.1.9) 

where v is the vertical projector induced by N and N„ is the Nijenhuis tensor 
associated to v, 

{X,Y) = [vX, vY] - V [vX, Y]-v [X, vY] + [X, Y], X,Y e X{E). 
In local coordinates, we set 

n 



-nfJdx' Adx^] 
2 



d 

dy"' 



and using (1.1.9), we obtain 

n{hX,hY) = -v[hX,hY], n{hX,vY) = n{vX,vY) = 0, (1.1.10) 
and then 

5N?' 

(1.1.11) 



^-^ 6x^ 6x^ 



The equality v = Id — h leads to = N/^, and we obtain 0, = — N/^. 

Proposition 1.1.6 The following equation holds 

[X^, Y^] = [X, Y]'' - n{X'',Y''). 

Using (1.1.10) follows a characterization of the integrability of the hori- 
zontal distribution. 

Theorem 1.1.1 The horizontal distribution of a nonlinear connection N is 
integrable if and only if the curvature vanishes. 

See [201 IH2] for the particular case of the tangent bundle. 
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1.2 Lie algebroids 



1.2.1 Cohomology 

Let M be a real, C°°-differentiable, n-dimensional manifold and (TM, ttm, M) 
its tangent bundle. 

Definition 1.2.1 A Lie algehroid over a manifold M is a triple {E, [•, ■]e, cr), 
where (E, vr, M) is a vector bundle of rank m over M, which satisfies the fol- 
lowing conditions: 

a) C°°{M) -module of sections T{E) is equipped with a Lie algebra structure 
[■,-]e. 

b) a : E ^ TM is a bundle map (called the anchor) which induces a Lie 
algebra homomorphism (also denoted a) from the Lie algebra of sections 
(T[E)^[-,-]e) to the Lie algebra of vector fields (x(-/Vf ), [•, •]) satisfying the 
Leibniz rule 

[Sl, fS2]E = f[si,S2]E + (tT(si)/)s2, Vsi, S2 G T{E), f G C°^(M). (1.2.1) 

Prom the above definition it results: 
1° {■,-]e is a M-bilinear operation, 
2° [■,-]e is skew-symmetric, i.e. 

[si,S2]e = -[S2, Si]e, Vsi, S2 G r(£'), 

3° {■■,-]e verifies the Jacobi identity 

[•Sl, [S2, SsIeIe + [•S2, [S3, Si]e]e + [^3, [si,S2\e]e = 0, (1.2.2) 

and fj being a Lie algebra homomorphism, means that 

cr[si,S2]E = [<y{si),cr{s2)]. (1.2.3) 

The existence of a Lie bracket on the space of sections of a Lie algebroid 
leads to a calculus on its sections analogous to the usual Cartan calculus on 
differential forms. In this paragraph we give only the relevant formulas for 
Lie algebroid cohomology we shall need later, and refer the reader to the 
monograph [73J for further details. 

If / is a function on M, then df{x) G -E* is given by {df{x), a) = cr{a)f, 
for Va G E^. For uj G /\'^{E*) the exterior derivative d^uj G /\''^^{E*) is 
given by the formula 
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k+1 

d^u}{si,...,Sk+i) =^ (-l)'+V(si)w(si,...,Sj,...,sjt+i)+ 

i=l 

+ ^ {-iy^^uj{[si^Sj]E,si,...,Si,...,Sj,...Sk+i). (1-2.4) 

l<i<j<k+l 



where Si € ^{E), i = l,k + l, and the hat over an argument means the 
absence of the argument. It results that 

{d^f = 0, (1.2.5) 

d^{ui A UJ2) = d^oji A 002 + (-l)'^^s'^ia;i A d^uj2, (1-2.6) 

The cohomology associated with d^ is called the Lie algebroid cohomology 
of E and is denoted by H*{E). Also, for ^ G r(£') one can define the Lie 
derivative with respect to ^ by 

C^ = i^od^ + d^ oi^, (1.2.7) 

where is the contraction with ^. 

1.2.2 Structure equations on Lie algebroids 

If we take the local coordinates (x*) on an open U C M, a local basis {s^} of 
the sections of the bundle 7r~^(C/) — )■ U generates local coordinates {x\y") 
on E. The local functions <T^(a;), L^^(x) on M given by 

9 



f^(sa) = f^a^i [sa,S/3]E = ^a^s^, i = l,n, a,/3,7 = l,m, (1.2.8) 

are called the structure functions of the Lie algebroid, and satisfy the struc- 
ture equations on Lie algebroids 

(a,/3,7) V / 

Locally, if / G C~(M) then d^ f = where {s"} is the dual basis of 

{sa} and if e G V{E*), 9 = ^^s" then 



d''0=U^^,-\e,Ll^]s-As^. (1.2.11) 



Particularly, we get 



1 



d^x' = ais"", d^s"" = -^L'^y A s^. (1.2.12) 
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Under a change of coordinates 



x*'=x*'(x*), i,i' = l,nonM, ^^ n ^o\ 

corresponding to a new base Sa = Sa', then the transformation rules of 
the structure functions are 



Some examples of Lie algebroids which will be used in this paper (see 
[75j for more examples and details) 

Example 1.2.1 The tangent bundle E = TM itself, with identity mapping 
as anchor. Then H*{E) = H*^ Rhami-^) Rham cohomology. With 

respect to the usual coordinates (x, x) , the structure functions are U-^, = 0, 
o"*- = (5*-, hut if we were to change to another basis for the vector fields, the 
structure functions would become nonzero. 

Example 1.2.2 Any integrable subbundle of TM is a Lie algebroid with 
the inclusion as anchor and the induced bracket. So, E = TT C TM, an 
involutive distribution associated with a regular foliation T , where one gets 
the tangential cohomology denoted H^^{M). 

Example 1.2.3 The cotangent bundle of a Poisson manifold E = T*M 
with (M, n) Poisson manifold. These carry a bracket characterized by the 
rule {df,dg} = d{f,g} and the anchor is the map 11* : T*M — )■ TM asso- 
ciated to the Poisson bivector field H. We obtain the Poisson cohomology 
denoted H^{M) fUU^ . 
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1.3 The prolongation of a Lie algebroid over the vector bun- 
dle projection 

Let {E, vr, M) be a vector bundle. For the projection vr : — )• M we can 
construct the prolongation of E (see [l6l [691 EOl EH ) • The associated vector 
bundle is {TE,tt2, E) where 

TE = U%nE, w e E, 

TwE = {{ux,Vw) G Ex X TyjE I cj{ux) = Tyjirivw), 7r{w) = x £ M}, 

and the projection tt2{ux,Vw) = t^e{'^w) = w-, where tte '■ TE E is the 
tangent projection. 

We also have the canonical projection vri : TE — t- E given by 7ri(n, = 
u. The projection onto the second factor : TE — )■ TE, a^{u,v) = v will 
be the anchor of a new Lie algebroid over the manifold E. An element of 
TE is said to be vertical if it is in the kernel of the projection tti. 

We will denote {VTE, ■K2\vte-' ^) vertical bundle of {TE, it2,E) and 

(t^Wte : VTE ^ VTE, 

is an isomorphism. If / € C°° (M) we will denote by and the complete 
and vertical lift to E of f defined by 

r{u)=a{u){f), f{u) = f{TT{u)), ueE. 

For s G T{E) we can consider the vertical lift of s given by 

for u G E, where 

is the canonical isomorphism. 

There exists an unique vector field s'^ on E, the complete lift of s satis- 
fying the two following conditions: 

i) is vr-projectable on a{s), 

ii) s'^{a) = LsOL, 

for ah a G r(^*), where a (n) = a(7r(n))(u), u^E (see [HI US]). 
Considering the prolongation TE of E over the projection vr, we may intro- 
duce the vertical lift and the complete lift of a section s £ ^(E) as 
the sections of TE — ^ E given by (see [70] ) 

s^{u) = (0, s^(n)), s^iu) = (s(7r(n)), s^(n)), u £ E. 

Another two canonical objects on TE are the Euler section C and the 
almost tangent structure {vertical endomorphism) J. 
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Definition 1.3.1 The Euler section C is the section of TE — t- E defined 
by 

C(n) = (0,0, Vu G E. 

Definition 1.3.2 The vertical endomorphism J is the section of the bundle 
(TE) e (TE)* E characterized by 

J(s^) = 0, J(s=) = s^, seT{E). 

The vertical endomorphism satisfies 

j2 = 0, ImJ = ker J = VTE, [C, J]te = -J- 

Definition 1.3.3 A section S ofTE E is called semispray (or second 
order differential equation -SODE) on E if 



J{S) = C. 

The local basis of T{TE) is given by {A'o, Va}, where [70 



d 

Xa{u) = ( Sa(vr(n)), cr'a-Q-i 



:i.3.i) 



(1.3.2) 



and {d / dx^ , d / dy°') is the local basis on TE. The structure functions of TE 
are given by the following formulas 



:i.3.3) 
:i.3.4) 



[Xa, Xis]rE = Ll^i^Xy, [Xa, Vis]rE = 0, [V^, V^Jte = 0. 
If y is a section of TE, then in terms of basis {^^a, Va} it is 

V = Z^'Xa + VVa, 
and the vector field cr^{V) G xi^^) has the expression 

The vertical lift of a section p = p'^Sa and the corresponding vector field are 
respectively 



d 
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The coordinate expressions of Euler section C and cr^(C) are 

and the local expression of J is given by 

j = ;f"®Va, (1.3.5) 

where {X" ,V'^} denotes the corresponding dual basis of {-^ajVa}. 
The Nijenhuis tensor of the vertical endomorphism 

Nj((z,'u;) = [Jz, Jw]e — J[Jz,w\e — J[z, Jw]e + J'^[z,w\e 

vanishes, and it results that J is a integrable structure. 

The expression of the complete lift of a section p = p'^Sa is 

p'^ = p-X^ + {ai^-L%(J')y'Va, 

and therefore 

^^p) = p^'^d^i + {^^^ - V) y g^- 

In particular 

si = Va, sl = X^- Li,fVp. 
The local expression of the differential of a function L on TE is 

jE t i dL dL 

and therefore, we have 

d^x' = aiX'^, d^y'^ = V^. 
The differential of sections of {TE)* is determined by 

d^ X'^ = -]^L%X^ A X^ , d^V'^ = 0. (1.3.6) 

In local coordinates a semispray has the expression 

S{x, y) = y'^X^ + 5"(x, y)V„. (1.3.7) 
The integral curves of o-^{S) satisfy the differential equations 

^ = 4(^K, ^=S^{x,y). (1.3.8) 
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If we have the relation 

[C,S]rE = S, 

then S is called a spray and the functions 5° are homogeneous functions of 
degree 2 in y". 

Under a change of coordinates = x* (x*), ?/" = A'^ (x*)y" on E the 
transformation rule of the coordinates on TE is given by 

X*' = x*'(x*), 
y°'=<(x%", 

and the rule of change of base is 

Va = A^Va'. 
The rule of change of dual base has the form 

X'' = A'^,X'^', 

dA°' 
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1.3.1 Ehresmann nonlinear connections on Lie algebroid TE 

As in the case E = TM I34[ HS] we can define the nonhnear connection. 



Definition 1.3.4 An Ehresmann nonlinear connection (or connection) on 
TE is an almost product structure N on 112 ■ TE E (i.e. a bundle 
morphism J\f : TE — )• TE, such that N"^ = Id) smooth on 7'-E'\{0} such 
that 

VTE = keT{Id + Af). 

If A/" is a connection on TE then HTE = kev{Id — J\f) is the horizontal 
subbundle associated to M and 

TE = VTE e HTE. 

Each p € T{TE) can be written as p = + where p^, p^ are sections in 
the horizontal and respective vertical subbundles. If p^ = 0, then p is called 
vertical and if p^ = 0, then p is called horizontal. A connection J\f on TE 
induces two projectors h, v : TE — )• TE such that h(p) = p^ and v(/j) = 
for every p G T(TE). We have 

h = l(/d + AA), v=l(/d-AA), (1.3.9) 

ker h = Imv = VTE, Imh = ker v = HTE. 
Locally, a connection can be expressed as 

M{Xa) = Xa-2M!^Vfs, M{Vp) = -Vp, (1.3.10) 

where J\fi = J\fa{x,y) are the local coefficients of A/". The sections 

6^ = h{X^) = Xa- M^Vp, (1.3.11) 

generate a basis of HTE. The frame {^q,, Vq} is a local basis of TE called 
adapted. The dual adapted basis is {A'",5V"} where 

Proposition 1.3.1 The Lie brackets of the adapted basis {da,Va} are 

(1.3.12) 

where 
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Proof. Using (1.2.1) and (1.3.11) we get 

/ dM^ ^dj\f^ -dMl dNl,\ 

If we insert Af^ = 6^ +N^Ve then the first relation from (1.3.12) is obtained. 
By direct computation the second relation is verified. □ 
We recall that the Nijenhuis tensor of an endomorphism A is given by 

Na(^, w) = [Az, Aw]te - A\Az, w]te - A[z, Aw]rE + A^[z, w]rE- 

Definition 1.3.5 The curvature of the connection M on TE is given by 
Q, = — Nh where h is the horizontal projector and is the Nijenhuis tensor 
ofh. 

Proposition 1.3.2 In local coordinates we have 

where T?.^^^ are given by (1.3.13) and represent the local coordinate functions 
of the curvature tensor in the frame /\^ TE* (g) TE induced by {A'q, Vq}. 

Proof. Since h^ = h we obtain 

U{z, w) = - [hz, hw]rE + H^z, w]rE + Hz, hwjrE - Hz, w]rE, 

Q{h.z,hw) = —v[h.z,hw]T-E, ^{h.z,vw) = Q{vz,vw) = 0, 
and in local coordinates we get 

which ends the proof. □ 
The curvature of the nonlinear connection is an obstruction to the in- 
tegrability of HTE, understanding that a vanishing curvature entails that 
horizontal sections are closed under the Lie algebroid bracket of TE. 

Remark 1.3.1 HTE is integrable if and only if the curvature ft = — Nh of 
the nonlinear connection vanishes. 

Let ^ a morphism of vector bundles E and E. We recall that the con- 
nections N on E and N on E are ^'-related if 

o AT = iVo ^f. 
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We consider the connections M on TE and N on TE which are cr^-related 
and a connection N on T^M which is o"*-related with N on TE and a - 
related with J\f on TE, where a : TE T^M is given by cr = cr* o u and 
cr* : — > T^M is the tangent apphcation of a. It follows 



Noa'^ = a'^ oM, N oa^ = a^oN,, Noa = aoM, (1.3.14) 

us 

by 



E „ TM „ 

Let us consider the adapted basis {Si, -^) of N and { 6 i, of N given 



f.= ^-Nf « 



and 



Therefore, we get 



TM _ d _ j^j d 



Theorem 1.3.1 The following relations hold 

MSi)=Si, |^+iV/=Ny^, (1.3.15) 

Proof. The first relation from (1.3.14) leads to the relation N((T^((5q,)) = 

(T^((5q) from which we get (T^{5a) = cr^ 6i and Afi = c^N^ . In the similar 
way the others relations are obtained. □ 

Proposition 1.3.3 For the curvature tensors of -related connections N 
and N we have the relation 

Kp = <4K^ (1-3-16) 

where 

is the curvature tensor of the nonlinear connection on TE. 
Proof. Using the relation = f a^f we obtain 

= {I m - i (n])) + NJ (^a^g - aj.^ j + L^^AT;, 

and from structure equations of the Lie algcbroid (1.2.10), the second term 
is NjcTgL^^ = —Me'L'^f^p, which concludes the proof. □ 
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Remark 1.3.2 A -related connection N on TE determines a connec- 
tion Af on TE with the coefficients Ma = cr^N^ and the curvature 7?.^^ = 

al^cTpE^Jj. The converse is not true because a is only injective. 

Let J be the vertical endomorphism. 

Remark 1.3.3 Let M be a bundle morphism of 1^2 '■ TE — >■ E, smooth on 
T-E\{0}. Then N is a connection on TE if and only if 

JM =J, MJ = -J. 

The proof proceeds as in the case E = TM and will be omitted. 

Definition 1.3.6 The torsion of a nonlinear connection N is the vector 
valued two form t = [J, h] where h is the horizontal projector and [■, ■] is the 
Frdlicher-Nijenhuis bracket. 

Definition 1.3.7 t is a semibasic vector-valued form. Its local expression 
is 

t=ltloX'^ AX^ <^V-y, (1.3.17) 



where 



Proof. We have 



[J, h]{z, w) = [Jz, hw]rE + [h-z, Jw]rE + J[z, w]rE - J[z, hi^jre - 
- J[hz,w]rE - Hz, Jw]rE - h[J2;, u;]7-E, 

and in local coordinates we get 

t{A:^,A:f^) =(^-^- LlA v„ t{x^,v^) = tiVa,v^) = o. 



□ 



Now, let us consider the linear mapping F : TE — >■ TE, defined by 

¥(hz) = -vz, ¥{vz) = hz, (1.3.19) 

for z G T(TE) and h, v the horizontal and vertical projectors of the nonlinear 
connection on TE. 

Proposition 1.3.4 The mapping F has the properties: 

i) ¥ is globally defined on TE, 

ii) Locally, it is given by 

W = -Va(^X'^ + Sa®6V'^, (1.3.20) 
Hi) F is an almost complex structure F o F = —Id. 
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Proof. It results by definition that F is globally defined and 

(F o F) (h^) = F (-vz) = -hz, (F o F) {\z) = F (h^) = -yz. 

In local coordinates we get F {5a) = —Va and F(Va) = 6a which ends the 
proof. □ 

Proposition 1.3.5 The almost complex structure is integrable if and only 
if the nonlinear connection is locally fiat, that is the curvature and torsion 
vanish. 



Proof. Let Np be the Nijenhuis tensor of the almost complex structure. 
We find 

Nf{5a,5p) = tl^5^-nl^V^, 



Nw{Sa, Vp) = -nl^S^ - tl^V^, (1.3.21) 

Nf{Va,V^) = --Nf{6a,6^). 

From (1.3.21) one reads immediately that Np = if and only if t = and 

n = o. □ 

A curve u : [to,ti] E is called admissible if a{u{t)) = c{t) where 
c{t) = 7r(ii(t)) is the base curve. A nonlinear connection on TE induces a 
covariant derivative of the sections defined locally as follows 



where p = p^Sa and rj = rj^'Sa- The derivative is linear in the first argument 
and it respects multiplication of second argument by real numbers, but 

not necessarily sum, except the case when the coefficients Ma arc the local 
coefficients of a linear connection. The linearity in the first argument permits 
us to define the derivative of a section rj G r(£') with respect to a e Ey, by 
setting 

where p G T{E) is satisfying p{u) = a. Also, the covariant derivative allows 
us to take the derivative of sections along curves. If we have a morphism of 
Lie algebroids ^ : F ^ E over the map ip : N ^ M and a section r] : N ^ E 
along (f, i.e r]{n) G -E^(n), n G iV, then rj can be written in the form 

p 

1=1 

for some sections {^i, ...,^p} of E and some functions {Fi, ...,Fp) G C°°{N) 
and the derivative of r] along ip is given by 

p 

1=1 
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where ap '^s the anchor map of the Lie algebroid F ^ N (see [30]). 

Let a : / — )• be an admissible curve and let 6 : / — )• be a curve in 
E, both of them projecting by vr onto the same curve 7 in M, 7r(a(t)) = 
7r(6(t)) = 7(i). Take the particular case of Lie algebroid structure TI — )• / 

and the morphism <I> : TI — )■ E, ^{t,t) =t 7(t) over 7 : / — )• M. Then one 
can define the derivative of b{t) along a{t) as Dd/dtbit)- In local coordinates, 
we obtain 

Definition 1.3.8 An admissible curve c{t) is a path (autoparallel) for non- 
linear connection N if and only if 

'^c{t)c{t) = 0. 

In local coordinates we get 

Prom the previous considerations we have: 

Proposition 1.3.6 An admissible curve c{t) in E is autoparallel for the 
nonlinear connection if and only if 



dx^ i Q, dy^ 



^ir, ^+AAiy" = 0, (1.3.22) 



where x' = x^{t) = x'{c{t)), = y°(t) = y"(c(t)), aj, = a'^{t) = cri{c{t)). 

Let be a nonlinear connection on TE, S' an arbitrary semispray on 
TE and h the horizontal projector of M. We consider S = hS' and for any 
other semispray S" on TE we have h{S' - S") = h((5'" - 5"")V„) = and 
it results that S does not depend on the choose of S'. We have 

JS = JhS' = JS' = C, 

so 5 is a semispray, which is called the associated semispray to M. 

Proposition 1.3.7 A nonlinear connection J\f and its associated semispray 
have the same paths. 

Proof. For the arbitrary semispray S' = y'^Xa + S'^^Va, the associated 
semispray of M is 

5 = h5' = y-^,-AA,VV/3, 

so 

= -M^y\ 

From (1.3.8) and (1.3.22) it results the conclusion. □ 



31 



Remark 1.3.4 If S is a semispray on TE, then we have 

J[S, Jz\te = -Jz, z £ T{TE). (1.3.23) 

Theorem 1.3.2 Let J be the vertical endomorphism on TE. If^ is a semis- 
pray then 

N = -jCsJ, (1.3.24) 

is a connection on TE. 
Proof. Since 

M{v) = -CsJ{v) = -[S, Jv\te + J[S, v\te 
using (1.3.23) we get 

JN{v) = -J[S, Jv\te + J'^[S,v\te = Jv, 

and 

j\fJ{v) = -[S, J^v]rE + J[S, Jv]rE = -Jv. 
By using the Remark 1.3.3 we get the proof of the theorem. □ 

Remark 1.3.5 The connection M = —CsJ is induced by the semispray S. 
Its local coefficients are given by 

Proof. By direct computation it results 

M{Xa) = -[S, J{Xa)]TE + J[S, Xo^TE 

and using (1.3.10) we obtain (1.3.25). □ 
Proposition 1.3.8 The torsion of the connection M = —CsJ vanishes. 
Proof. We have 

t = [J, h]TE = \ ([J, Id\TE + [J, -[e, J]te\te) = ^[J, [J, ^]te]te- 

Using Jacobi identity wc obtain tliat t = 0. Also, if we use (1.3.25) into 
(1.3.18), by direct computation, the same result is obtained. □ 
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Proposition 1.3.9 The associated semispray of N = —CsJ is given by 

^{S-[S,C]te). 
Proof. The associated semispray is 

hS = ^S + ^M{S) = \{S- [S, JS]te + J[S,S]te) = \iS- [S,C]te). 

□ 
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1.3.2 Lagrangian formalism on Lie algebroids 

We consider the Cartan 1-section 

which, in local coordinates is 



e, = — (1.3.26) 



The differential 9l is the Cartan 2-section 

and from the local coordinate expression of 0^ we get 



But 

and it results [70] 



(1.3.27) 

The function L is said to be a regular Lagrangian if is regular at every 
point as a bilinear form. Let us consider the energy function given by 

dej dL 

and the symplectic equation 

isUL = -d^EL, S e T{TE). (1.3.28) 
In local coordinates, considering the section 

we obtain the equations 
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The equality of the V" components yields 

5a/3(/-/)=0, 

and the regularity of the Lagrangian implies = which means that 
S = y'^Xa + S"Va is a semispray. Analogously, the equality of the Af" 
components leads to the equation 

d^L dL dL 

and the regularity condition of the Lagrangian determines the components 
of the semispray 



'dx^ ""dx^dyl^^ '^"^ dy^ 



where Qapg^^ = Sa- From (1.3.29) and (1.3.25) it results: 

Corollary 1.3.1 For a regular Lagrangian L, there exists a nonlinear con- 
nection N with the coefficients given by 



1 / dS^ 

2 \dr 

where ^ 

and will be called the canonical nonlinear connection induced by a reg- 
ular Lagrangian L. 

If (x*) are coordinates on M, {s^} is a local basis of T{E), {x'\y°^) are 
the corresponding coordinates on E and ^{t) = (x*(t), ?/"(t)) then, 7 is a 
solution of the Euler-Lagrange equations if and only if [125] 

dfdL\ i dL je ^fidL 
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1.3.3 Homogeneous connections 
Definition 1.3.9 The morphism 

T=\ccM, 
is called the tension of the nonlinear connection. 
In local coordinates we get 

T{Xa) = (m^ - l^y^j Vp, r(V.) = 0, 

and it results 

T=i^^-^-^f'^X'^®Vp. (1.3.32) 

It is obvious that T is vanishing, if and only if the nonlinear connection is 
homogeneous of degree 1 with respect to y". 

Proposition 1.3.10 If S is a spray then N = —CsJ is a homogeneous 
nonlinear connection. 

Proof. Using (1.3.25) we get 

' ~ \ dy- +^ dy-dyP) ^^^^ 
But 5 is a spray and it results that is homogeneous of degree 2, there is 

25- - 

and 



Qya dy'^dyl^' 

therefore, the tension vanishes. □ 
Definition 1.3.10 The strong torsion T of M is given by 

T = i^t- T, 

where T is the tension, t is the torsion of M, and i^ is the contraction with 
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Locally, we obtain 

'dN7 

dy 



Proposition 1.3.11 The strong torsion T of a nonlinear J\f connection 
vanishes if and only if the torsion and the tension of N vanish. 

Proof. If T = then we have 

ATT = I y/3_J^7 



and it results 



which yields t((5a, 5^) = and T = 0. □ 

Definition 1.3.11 A function :£'—)• [0, oo] which satisfies the following 
properties 

1) isC°° on E\{0} 

2) T{Xu) = XT{u) for X> andue E^. X e M. 

3) For each y G Ex\{^} the quadratic form 

Qafi\x,y) 



2 dy'^dyf^ ' 

is positive definite, will be called the Finsler function on a Lie algebroid. 

If we insert C = ^J^"^ = ^gajsy'^y^ into the expression of semispray 
(1.3.29) we obtain 



Corollary 1.3.2 A homogeneous nonlinear connections has the coefficients 
given by 

with 

cS ^ SB f i ^9/37 , i ^9ali i dga-y , r£ , re re A a T 

^ =2^ y^dS' + " + ^^"^^^ + " ^^'^^'^ J ^ y ■ 

(1.3.33) 

and is called the canonical nonlinear connection associated to a Finsler func- 
tion. 
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Remark 1.3.6 In the particular case of the standard Lie algebroid E = TM 
and a = Id the Cartan nonlinear connection is obtained. 

We consider the canonical nonlinear connection and ||y||^ = Qapv'^y^ = 
is the square of the norm of the Euler section. The almost complex 
structure characterized by (1.3.20) does not preserve the property of homo- 
geneity of the sections. Indeed, it applies the 1-homogeneous section 5a onto 
the 0-homogeneous section Vq, a G l,m. We define a new almost complex 
structure Fq : TE — >■ TE given by 

M^a) = --Va, Fo(Va) = |;5a, a > 0. 

It is not difficult to prove that Fq = —Id and Fq preserves the property of 
the homogeneity of the sections. 

Theorem 1.3.3 The almost complex structure Fq is integrable if and only 
if the following relations hold 

(1.3.34) 

where ya = g^^y^ , a,P,-f = l,m. 

Proof. For the Nijenjuis tensor Nutq we have 

- {tip + 2^ {H^)^c. - 6a{F^)5})) V^, 

NFo(<5a,<^/3) = -f^NF„(Va,V;3)- 

It follows that Nfq = if and only if the relations (1.3.34) are satisfied. □ 
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1.3.4 Linear connections on Lie algebroids 

A linear connection on a Lie algebroid {E, [,]e,<^) is a map 

V : r{E) X r{E) r{E), 

which satisfies the rules 

i) Vp+^r] = VpTj + V^rj, 

ii) Pp(?? + uj) = VpT] + VpOj 

iii) VfpTi = fV^^, 

iv) Vp{fn) = ia{p)f)v + fV^r, 

for any function / G C°°(M) and p,r],oj e T{E). 

For p,ri & r(-^) the section VpT] G ^(E) is called the covariant deriva- 
tive of the section rj with respect to the section p. Let A/" be a nonlinear 
connection. We have: 

Definition 1.3.12 A linear connection V on Lie algebroids is called N— 

linear connection if 

i) D preserves by parallelism the horizontal distribution HTE. 

ii) The tangent structure J is absolute parallel with T>, that is T>J = 0. 

Consequently, the following properties hold: 

{Vpri^y = 0, {Vprff = 0, Vph = 0, VpV = 0, 
Vp{Jrl-) = JiVpn"^), Vp{Jrf) = J{Vprf). 

If we denote 

2^> = 25phr?, V;ri = Vp.ri, 
then the following decomposition is obtained 

Vp = v^p + v;, peViE). 

We remark that and are not covariant derivative, because 'Dpf = 
(^{P^)f (^{P)f, T^pf = (^{P^)f ^ (^{p)f, but, it stiU preserves many proper- 
ties of v. Indeed, and T>^ satisfy the Leibniz rule, and and P"^ will be 

called the h— covariant derivation and v— covariant derivation, respectively. 
Using the fact that T>s„ = , I?Vq = T^Va ^et 

Proposition 1.3.12 In the adapted basis {6a, Va} a M— linear connection 

can be uniquely represented, in the form 
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The system of functions {F2p{x,y), C2^(x, y)) represent the local coeffi- 
cients of h— covariant derivation and of v— covariant derivation, respectively. 
Under a change of coordinates (1.2.13) the coefficients satisfy the rules 

_ (-11' AOi' aP' A-y 
Let us consider a c?— tensor T in the local adapted basis, given by 

A d— tensor means a tensor on TE, whose components, under a change 

of coordinates on TE behave like the components of a tensor field on the 
base manifold E. Its covariant derivative with respect to ^ = + = 
^"Sa + C^Vp is given by 

where we have the h— covariant derivative T>^T = ^'^T^^^^Sa(S>Vi3^X'^ CSiSV^ 
with 

and is the operator of h— covariant derivative. 

The v-covariant derivative of T is VJT = ^'^T^^^/Ja V/j Af ^ SV^ with 



7£ IK Qyi^ TK 'ye ' ^rn 'ye ^ek 7T re 



where is the operator of v— covariant derivative. 
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1.3.5 Torsion and curvature of a J\f— linear connection 
The torsion tensor of a A/"— linear connection is defined as usual 

As in the case of tangent bundle we have: 

Proposition 1.3.13 The torsion of a M— linear connection is completely 
determined by the following five components 

' hT(he, he.) = V\\iu: - Pjjhe - h[h^, ^\rE. 

vT(h^,ha;) = -v[h^, ha;]r£;, 
< hT(hC,va;) = -P>e-h[he,vu;]ri=;, (1.3.35) 

vr(h^, vw) = V\\u - v[h^, vwJte, 
. vT(v^, vw) = P^vo; - P^vC - vK, va;]rE. 

With respect to the adapted basis the components of torsion are given 

by 

< hT{Vp,5a) = Cl^5^, 

^ vTiVf,, Va) = Sl^V, = (C^^ - CjJ) V,. 

The curvature of a A/"— linear connection is defined by 

Proposition 1.3.14 The tensor of curvature has three essential compo- 
nents 

R{5j,Si3)5a = R^^^Ss, 
< R{V^,5^)6a = P^p^S„ (1.3.36) 
[ i?(V„V^)(5„ = 5^^^4, 

given by 

I 171(5 pS TTip jP 171(5 \ -1->£ 
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dKin ■ dC^s: dC^s: dN% , 



Proposition 1.3.15 The Ricci identities have the following form 

-\'Ci \ra vP t^P vct n-}p v-a i 

p dxi \ T dx^ ^ dx^ i dx^ ' 



X°'/p/^ - X"" 1^1 = Sp/^^XP - SPp„X"-/p. 

Definition 1.3.13 A M— linear connection is called of Cartan type if 

V\C = 0, VJC = v^, (1.3.37) 

where C = y"Va is the Euler section. 

By direct computation, it results that a A/^— linear connection D is of 
Cartan type if and only if 

A^^ = i^e>^ y'Q"^ = 0. (1.3.38) 

Introducing these relations into the coefficients expression of the curvature, 
we obtain the following result: 

Proposition 1.3.16 A J\f— linear connection of Cartan type has the prop- 
erties 

TZp^ = R^l^^y^ - L'^^, Pp^ = P^p^y^, S^^ = S"p^y^. (1.3.39) 
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1.4 Dynamical covariant derivative and metric non-linear 
connection on Lie algebroid TE 

In this section we will introduce the notion of dynamical covariant derivative 
on Lie algebroids as a tensor derivation and study the compatibility between 
nonhnear connection and a pseudo-Riemannian metric. 

Definition 1.4.1 A map V : %{TE\{{)}) %{TE\{{)}) is said to be a 
tensor derivation on 7'-E'\{0} if the following conditions are satisfied: 
i) V is M-linear 

a) V is type preserving, i.e. V(T^(T-E\{0}) C 'Zl{TE\{0}), for each 
{r,s) G N X N. 

Hi) V obeys the Leibnitz rule V(P(X'S') = VP(X'S' + P(X' VS, for any tensors 
P,S onTE\{0}. 

iv) V commutes with any contractions, where T*(T-E'\{0}) is the space of 
tensors on TE\{Q}. 

For a semispray S we consider the M-hnear map 
V^:T{TE\m ^T{TE\m, 

given by 

Vop = h[S, hp]rE + v[5, vp]rE, Vp G TirE\{0}). (1.4.1) 
It results that 

Vo(/p) = S{f)p + fVoP, V/ GC^iE), p€ r{TE\{0}). (1.4.2) 

Any tensor derivation on 7'-E\{0} is completely determined by its actions on 
smooth functions and sections on 7'£'\{0} (see |117j generalized Willmore's 
theorem, p. 1217). Therefore there exists a unique tensor derivation V on 
TE\{0} such that 

V \c"°(E)= <S, V |r(r£;\{o})= ^o- 

We will call the tensor derivation V, the dynamical covariant derivative 
induced by the semispray S and a nonlinear connection A^. 

Proposition 1.4.1 The following formulas hold 

[S,Vp]TE = -5p - [M^ + -g^j V„, (1.4.3) 
[5, 5^]rE = {M^ - L%f) 6a + n}V^, (1.4.4) 



where 

n = -^^§1 - +M^M2+Mf^ +M]Ll,y-. (1.4.5) 
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The action of the dynamical covariant derivative on the Berwald basis 
is given by 

V6p = h[S, 5p]TE = {M^ - L%f) 60,. 

It is not difficuh to extend the action of V to the algebra of tensors by 
requiring for V to preserve the tensor product. For a pseudo-Riemannian 
metric g on VTE (i.e. a (2, 0)-type symmetric tensor g = ga[s{x, y)V" (8) 
of rank m on VTE) we have 

(V5)(/>i,/>2) = S{g{pi,p2)) - g{Vpi,p2) - g{pi,Vp2), (1.4.6) 
and in local coordinates we get 

(1.4.7) 

Definition 1.4.2 The nonlinear connection N is called metric or compati- 
ble with the metric tensor g if \/g = 0, that is 

S{g{pi,p2)) = g{Vpi,P2) + 9{P1,VP2). (1.4.8) 

If 5 be a semispray, M a nonlinear connection and V the dynamical 
covariant derivative induced by (5,AA), then we set: 

Proposition 1.4.2 The nonlinear connection N with the coefficients given 
by 

M^=M^-\g''^g,p,, (1.4.9) 
is a metric nonlinear connection. 

Proof. Since A'^ are the coefficients of a nonlinear connection and 

9°''^ 9'yl3/i are the components of a tensor of type (1,1) it results that N"^ are 
also the coefficients of a nonlinear connection. We consider the dynamical 
covariant derivative induced by (5, N) and we have 

{Vg){Vo.,Vp) = S{g^p)+g,p(^Kf2 + ^^+g,c.{KfJ + ^ 

= S{g^p) + g,p {m2 + IP) + g,a [mJ + ^) - 

~9-if}\9^''9ea/-9^J^9^''9ef}/ = 

1 1 

= 9ali/ - 29al5/ - 29al3/ = 0, 

that is the connection M is metric. □ 
For the particular case of the tangent bundle see [191 EI] 
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1.4.1 The case of SODE connection 

A semispray {SODE) given hy S = y°'Xa+ S'^Va determines an associated 
nonlinear connection 

M = -CsJ, 



with local coefficients (1.3.25) 

Proposition 1.4.3 The following equations hold 

[S, Vp]rE = -6p + [M^ - L%f) Va, (1.4.10) 

[S, 5p\rE = {Mj^ - L%y') Sa + 7^^V«, (1.4.11) 

where 

n = -^"P^ - S^^p) - K-^J + (LJi^K + L^-^M- (1-4-12) 

The dynamical covariant derivative induced by S and associated nonlin- 
ear connection is characterized by 

VVf, = v[5, V^]rE = {M^ - L%y^) = - 2 + L%y') (1-4.13) 

V5p = h[5, 5p]TE = [N^ - L%f) 6^. 

9aP/ ■= (V5)(Va, V/s) =S{galj) - g^pM2 - g^aMJ + {g^H^Le + giaL'le) 

(1.4.14) 

which is equivalent to 

(1.4.15) 

Definition 1.4.3 The Jacobi endomorphism is given by 

$ = v[5, hp]rE- 

Locally, from (1.4.11) we obtain that $ = 7^^V„0A''', where IS given 
by (1.4.12) and represent the local coefficients of the Jacobi endomorphism. 

Proposition 1.4.4 The following result holds 

$ = is^ + v[v(S, hplrE- 
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Proof. Indeed, $(^) = v[5,hp]7-£; = v[h5, h/o]7-£; + v[v5, hpjT-^ and 
n{S, p) = v[h5, \ip\rE, which yields $(p) = Vl{S, p) + v[vcS, \ip\rE- □ 

If 5 is a spray, then the coefficients 5" are 2-homogeneous with respect 
to the variables and it results 

flca 

5 = h5 = y«5„, v<S = 0, M^ = ^y' + L%f, 
which yields 

$ = ^5^^, (1.4.16) 

and locally we get 

ni = 7^°^y^ (1.4.17) 

which represents the local relation between the Jacobi endomorphism and 
the curvature of the nonlinear connection. 
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1.4.2 Lagrangian case 

Let us consider a regular Lagrangian L on E, that is the matrix 



9al3 



dy^dyP ' 



has constant rank m. The symplectic structure induced by the regular 
Lagrangian is [70] 



Let us consider the energy function given by 
then the symplectic equation 

and the regularity condition of the Lagrangian determine the components 
of the semispray (1.3.29) 

where gajsg'^^ = SZ. 

The connection J\f determined by this semispray is the canonical nonlinear 
connection induced by a regular Lagrangian L. Its coefficients are given by 



I 

' dx^dy'^ dx^dy/^ dy' 

(1.4.18) 



Theorem 1.4.1 The canonical nonlinear connection N induced by a regu- 
lar Lagrangian L is a metric nonlinear connection. 



47 



Proof. Introducing the expression of the semispray (1.3.29) into the 
equation (1.4.15) we obtain 

1 / dg'" dg^^\ f i dL ^ d'^L g . dL 



Ua' V^^a^ ) 



By direct computation, using the equaUties 

it results (V5')(Va, V/j) = 0, which ends the proof. □ 

Theorem 1.4.2 The canonical nonlinear connection induced by a regular 
Lagrangian is a unique connection which is metric and compatible with the 
symplectic structure ojl, that is 

V5 = 0, (1.4.19) 

ul{^p, hi/) = 0, Vp, u G r(r£;\{o}) (1.4.20) 

Proof. Using the equation V" = (5V" — M^X^ it results 

= «(*v''-<.v^)A.v«+i (^< - ^4 - 0l;,) .v*a^« 
= w*v*A^«+l (^,, - Af,. + - - 0i«,) A-^A^", 

where A/'a^ := ga-y^J- We have that a;L(hp, hi/) = if and only if the second 
part of the above relation vanishes, that is 
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It result that the skew symmetric part of Nap is uniquely determined by the 
condition (1.4.20). The symmetric part of Nap is completely determined by 
the metric condition (1.4.19). Indeed 

S{gap) = g-tpN2 + g^aNJ - {g^pL1^^ + g^aL^^^y" 

= Npa + NaP - {g^l3L''L + 9-taLl^ y'' 

that is 

2AA(„^) = Sigap) + {g-ipLl^ + g-taL]^e) y'- 

The equations (1.4.19) and (1.4.20) uniquely determine the coefficients of 
the nonlinear connection 

N; = g^"Nap=g^"{N(aP)+N[aP]) 



2^ 

+^5^" (g-rpLl^ + 9jaLl^) y' 

Conversely, introducing (1.3.29) into (1.3.25) we have (1.4.18) which ends 
the proof. 

Remark 1.4.1 The invariant form of Helmholtz conditions on Lie alge- 
broids is given by: 

Vlg{v,e)=Vlg{v,p), 

V<7 = 0, (1.4.21) 

g{^p,v) =g{^i^,p), 

for v^p^O S r(£'), where is v-covariant derivative. 
In local coordinates we obtain the following equations 

dgaP _ dgae 

S{gap) - g-ypNl - g^aNj = y" {g^pLJ^ + g^aLl^^ 
5-7 [^'p% + SNJ + NINJ - {L%NJ + LlNf,)f^ = (1.4.22) 

9P, {^<% + SNl + N^N] - {Li^N] + LlNi)f 

In the case of standard Lie algebroid (TM, [■, •], id) we obtain the classical 
Helmholtz conditions |113j . 
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1.5 The prolongation of a Lie algebroid over the vector bun- 
dle projection of the dual bundle 

Let T : E* ^ M he the dual bundle of tt : — )• M and {E, [•, •]_£;, cr) a Lie 
algebroid structure over M. One can construct a Lie algebroid structure 
over E*, by taking the prolongation of {E, [•, •]e,o') over t : E* ^ M (see 
[36] . [M] . [50]). This structure is given by: 

• The associated vector bundle is [TE* ,ti, E*) where 

TE* = UTu-'E*, u* £ E* 

with 

Tu*E* = {{u^,Vu*) G E^ X Tu'E*\a{u^) = Tu*tM,t{u*) = x E M}, 

and the projection ti : TE* — )• E'*, Ti(Ma;,f„*) = n*. 

• The Lie algebra structure [•, ■]te* on TiJ'E*) is defined in the following 
way: if pi,p2 G TiJ'E*) are such that Pi{u*) = {Xi{T {u*)),Ui{u*)) where 
Xi G r{E), Ui G x{E*) and ^(^^(r (n*)) = T^^riUi (u*)), i = 1,2, then 

[pi,P2]ri^*(^^*) = i[X,,X2]TE*ir iu*)),[UuU2]TE'{n*)). 

• The anchor is the projection : TE* — t- TE*, a^{u,v) = v. 
Notice that if Tr : TE* E,T t{u,v) = u then (VTE* ,Ti\vrE* , E*) 

with VTE* = KerTr is a subbundle of (TE* ,ti,E*), called the vertical 
subbundle. If {q^,fia) are local coordinates on E* at u* and {sa} is a local 
basis of sections of vr : — )• M then a local basis of r(7'i^'*) is {Qa,^"} 
where 

QaK) = (^a(r(^*)),<^|n.) , ^'"(n*) = (o,^ln') . (1.5.1) 
The structure functions on TE* are given by the following formulas 

[eo.esiTE. "ij^e,, [s-,?'"]™. = o, [7:'",p'']rE. = o, (1.5.3) 

and therefore 

d'^q' = oiQ', d^Pa = Va, (1.5.4) 

where {Q"',Va} is the dual basis of {Q«,P"}. Also, if p = p°Q^ + p^T'" is 
a section of T-E* , then 

1/ X i a d , d 

^ {P) = <P TT- + Pa 
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If u* G E* and (n^,Uu*) € x Tu*E* then 

is called the Liouville section. The canonical symplectic section coe defined 

by 

OJE = -d^0E, 

and it results that is a nondegenerate 2-section and d^ujE = 0. 

On E* we have the similar concept of the vertical lift to that in ii^. If 
a G T(E*) we can define the vector field on E* as follows 

a''{u*) = a{T{u*))l., u* G E* , 

where 

is the canonical isomorphism between the vector spaces and T„* {E*^^^-^) 

Also, if X is a section of vr : — )• M, there exists a unique vector field X*'^ 
on i?*, called the complete lift of X to E* , satisfying the two following 
conditions: 

(i) X*^ is r— projectable on cr(X), 

(ii) = C^, for all y G r(S) (see [H]). 

A A 

If X is a section of E then X is the linear function XG C^{E*) given by 

X{u*) = u*{X{t{u*))), 

for all u* £ E* . Now, we may introduce the vertical lift and the complete 
lift X*^ of a section a G T{E*) and a section X G r(i?) as the sections of 
TE* given by 

a^(n*) = (0,a^(n*)), X*'' = {X{t{u*), X*%u*)), u* £ E* . 
The other canonical object on TE* is Liouville- Hamilton section C given by 

C(n*) = (0,n™), u*£E*. 
In local coordinates it follows that the Liouville section is given by 

eE = l^aQ'', (1.5.5) 

and we obtain 

U;E = Q"AVa + \lJ^aL%Qf' A Q\ (1.5.6) 

The Liouville- Hamilton section C has local expression 

C = HaV. (1.5.7) 
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If ^ is a section of E*, 6 = 9as", and X is a section of E, X = X°'sa then 
the vertical and complete lifts have the expression 

The Liouville-Hamilton section on TE* measures the homogeneity of the 
functions and sections. A function / G C°°{E*) is said to be homogeneous 
of degree r G Z if 

Ccf = rf, 

where Cc is the Lie derivation with respect to the Liouville-Hamilton section 
on the Lie algebroid. A section p of TE* is said to be homogeneous of degree 
r G Z if 

jCcP = rp. 

We remark that VTE* is Lagrangian for we, i.e. oje{pi,P2) = 0, for every 
vertical sections p\,p2 € TiyTE*). 
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1.5.1 Ehresmann connections on the Lie algebroid TE* 

Definition 1.5.1 The Ehresmann nonlinear connection on TE* is an al- 
most product structure N on ti : TE* E* (i.e. a bundle morphism 
N : TE* TE*, such that = Id) smooth on TE*\{Q] such that 
VTE* = keT{Id + M). 

If is a connection on TE* then HTE* = ker(Id—J\f) is the horizontal 
distribution associated to M and 

TE* = VTE* ®HTE*. 

Each p G T{TE*) can be written as p = + p" where p^, p" are sections 
in the horizontal and respective, vertical subbundles. If p'^ = then p is 
called vertical and if = then p is called horizontal. A connection M 
on E* induces two projectors h, v : TE* TE* such that h{p) = p^ and 
v{p) = p" for every p G V{TE*). We have 

h=^-{Id + N), v = ^{Id-N), 

ker h = Imv = VTE* , Imh = kerv = HTE* , 
h^ = h, = V, hv = vh = 0, h + v = Id. 
Locally, a connection can be expressed as 

^^{Qa) = Qa + v^a^r^ H{r'') = -v', 

where J\fa/3 = -f^af^iQ, 1^) are the local coefficients of A/". The local coordinate 
expression of M is 

The local sections P", (a = l,m) define a local frame of VTE*, and the 
sections 

S*a = {Qa)'' = Qa+Ma0r^ (1-5.8) 

generate a local frame of HTE*. The frame {(5* , 7^"} is a local basis of TE* 
called adapted to the direct sum decomposition. The dual adapted basis is 
{Q",SVa} where 

SVa=Va-Maf}Q^. (1.5.9) 

It results that in the adapted basis the expression of Ehresmann connection 
becomes 
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Definition 1.5.2 A connection M is called symmetric if HTE* is Lagrangian 
for UJE- 

By a straightforward computation, using (1.5.8) we get 
and it results that M is symmetric if and only if 

^fal3 - ^Pa = ^7-^I;3- (1.5.10) 

Proposition 1.5.1 With respect to a symmetric nonlinear connection, the 
canonical symplectic structure cve can be written in the following form 

UE = Q" ASVa + f^aL'^^Q" /\Q^- (1.5.11) 

Proof. Using (1.5.6) and (1.5.9) we get 

which ends the proof. □ 
Proposition 1.5.2 The Lie brackets of the adapted basis {(5*,7^"} are fSUjj 

Proof. Using (1.5.8) we obtain 

\K. h\TE^ - \^^^ - + ^-^^-^^^^ ) ^ ' + 

and putting Qe = 6*, - MejV^' we get [5* , S^jrE* = Ll^S; + TZ^i^^V^' . □ 
The curvature of a connection J\f on TE* is given by = — N/^, where 
h is the horizontal projector and N/^ is the Nijenhuis tensor of h, given by 

Nh{9,p) = [he,hp]rE' - h[he,p\TE' - h[e,hp]rE' + h^[9,p]rE'. 

Proposition 1.5.3 In local coordinates we get 



^^a/37Q" AQ^^T'^, (1.5.13) 



0, — —-/-Ca/S-y 

where TZa/S'y is given by (1.5.12) and are called the coefficients of the curva- 
ture tensor of M . 
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Proof. Since = h we obtain 

n{hpi,hp2) = -v[hpi,hp2]TE*, 0,{hpi,vp2) = 0,{vpi,vp2) =0, 
and in local coordinates we get 

n{si,s*^) = -v[si,s}]rE* = -n^0jV\ 

which concludes the proof. □ 
Proposition 1.5.4 The curvature satisfies the Bianchi identity 

Proof. By direct computation, using the relation (1.5.12) and structure 
equations given by (1.2.9), (1.2.10). 

The curvature is an obstruction to the integrability of HTE*, under- 
standing that a vanishing curvature entails that horizontal sections are 
closed under the Lie algebroid bracket of TE*. We have: 

Remark 1.5.1 HTE* is integrable if and only if the curvature vanishes. 

Also, the integrability conditions for the almost product structure is 
given by the vanishing of the associated Nijenhuis tensor N/^-By a straight- 
forward computation we obtain 

N^(P",P^) = 0, Nj^i5*^,V^) = 0, Nj^{S*^,6}) = ma0^V\ 

Thus 

Nat = -27^«^^Af" AXl^ ^ 

and it results that the distribution HTE* is integrable if and only if the 
almost product structure M is integrable. 

We consider the connections M on TE* and N on TE* which are a^- 
rclatcd and the adapted basis {6i, ■^-) of N given by 5i = ^ + Ni^-^;-, 
where N^p are the coefficients of A''. 

Theorem 1.5.1 The following relations hold 

(^HS*a)=<^u Ma0 = C7iNif3, TZafS^ = a'y^Rij^, (1.5.14) 

where 

is the curvature of the nonlinear connection N on TE* . 
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Proof. Since J\f{a'^{6l,)) = ct^(5*) the relation N o = oA^ leads to 
N{(7'^{5*J) = (T^{6*J. But iV(c7i((5*)) = 2ai6i-a^{S*^) which concludes that 
C7^((5*) = al^Si, from which we easily obtain Mafi = cra^ifi- By straightfor- 
ward computation we get 

and using (1.2.10), the second term is given by Nj^a^L^^^ = — iVg-yL^^, which 
ends the proof. □ 

Remark 1.5.2 A a^-related connection N on TE* determines a connection 
N on TE* with the coefficients given by 

and curvature 

Conversely, it is not true, because a is only injective. 
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1.5.2 Regular sections and connections 

Definition 1.5.3 An almost tangent structure J on TE* is a bundle mor- 
phism J : TE* TE* of n : TE* E* of rank m, such that = 0. An 
almost tangent structure J on TE* is called adapted if 

ImJ = keTj = VTE*. 

Locally, an adapted almost tangent structure is given by 

J = tapQ''®V^, (1.5.15) 

where the tensor tajsix, fi) is nondegenerate. 

Now, we find the integrability conditions for tlie adapted almost tangent 
structure. From 1501 it results 



Proposition 1.5.5 J is an integrable structure if and only if 



(1.5.16) 



where t"-'^t^p = 5'^. 

Proof. J is integrable if and only if the associated Nijenhuis tensor 

Nj(P,^) = [Jp,Jy]TE* - J[Jp,y]TE*-J[p,Jy]TE', 
vanishes. This is locally equivalent with the relations 

Therefore is integrable if and only if 

dt^ _ ^ dtae 

that is equivalent to (1.5.16). 



ta-y n tRy 

Ofiy Oily 



Definition 1.5.4 An adapted almost tangent structure J on TE* is called 
symmetric if 

oje{Jpi,P2)=^e{JP2,Pi)- (1.5.17) 
Locally, this requires the symmetry of the tensor tap- 
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Remark 1.5.3 If g is a pseudo-Riemannian metric on the vertical bundle 
VTE* then there exists a unique symmetric adapted almost tangent structure 
on TE* such that 

g{Jp, Jv) = -uje{Jp, v), (1.5.18) 
and we say that J is induced by the metric g. 

Locally, if 

g{q,fi) =g"^Va(^Vfs, 
then the relation (1.5.18) implies t°'^ = g'^^ . 

Remark 1.5.4 Any symmetric adapted almost tangent structure J on TE* 
induces a pseudo-Riemannian metric on the vertical bundle VTE* as defined 
by (1.5.18). 

Definition 1.5.5 The torsion of a connection M is the vector valued two 
form T = [J', h] , where h is the horizontal projector and [J', h] is the 
Prdlicher-Nijenhuis bracket 

[J, h]{X, Y) = [JX, hY]rE* + [hX, JY]te* + J[X, Y]rE* - J[X, hY]rE* 
-J[hX,Y]rE^ - h[X,JY]rE* - h[JX,Y]TE*, 

The torsion T is a semibasic vector-valued form. Its local expression is 



TaM = tae r^,, ^ - ^/3e r,,"^ + Saitp-y) - S*p{taj) - LaffU-y. (1.5.19) 



where 

Next, let us consider the linear mapping T : TE* — )■ TE* given by 

J^{hp) = Jp, F{Jp) = -hp, (1.5.20) 

where p G r(T-E*), and h is the horizontal projector induced by the nonlin- 
ear connection. 

Proposition 1.5.6 The mapping T has the properties: 

a) J- is an almost complex structure J- o T = —Id. 

b) Locally it is given by 

J^ = tapV^®Q"-t"^6*^^5Vp. (1.5.21) 
Proof. It results by definition that 

o F){hp) = HJ{P)) = -hp, {T o T){Jp) = H-hp) = -Jp, 

T{dl)=tafiP^ J-(P°) = -rt^, 

which concludes the proof. □ 
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Proposition 1.5.7 The almost complex structure is integrable if and only 
if the torsion and curvature of the connection satisfy the equations 

Tq,^^ = 0, 7Cq,/3^ = tae-;\ — • (1.5.22) 

Proof. Let Njr be the Nijenjuis tensor of the almost complex structure. 
We get 

atT , ■ , g-'-y -« -J- dt jS'y ^"^(yy 

= Ta/Set ^al3{-y) = ^"^~qJ^ ~ ^^^^^^ ~ 



where 



which ends the proof. □ 

Remark 1.5.5 Let N he a bundle morphism of ti : TE* — )• E* , smooth on 
T-E'*\{0}. Then M is a connection on TE* if and only if there exists an 
adapted almost tangent structure J on TE* such that 

JM = J, MJ = -J. (1.5.23) 

Definition 1.5.6 Let J be an adapted tangent structure on TE* . A section 
p ofTE* is called J'— regular if 

J[p,Ji^]te* = -Ji^, (1.5.24) 

for every section v of TE* . 

Locally, the section p = (f^Qa + PjiT^ is J"— regular if and only if 



t 



dpa 



where f^^t^^ = 5^. 

We have to remark that if the equation (1.5.24) is satisfied for any section 
G T{TE*) with ronA;[t"^] = m, then J is an integrable structure. Indeed, 
we have 

dtf^ _ d^^P _ Q^e^ _ dt^P 

and using (1.5.16) it follows that J' is integrable. 
From ^50j we have: 
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Theorem 1.5.2 Let J he an adapted tangent structure on TE*. If p is a 
J— regular section ofTE* then 

N=-CpJ, (1.5.25) 

is a connection on TE* . 
Proof. Since 

M{v) = -CpJ{v) = -[p,Jv]rE'+J[p,v]rE*, 

then 

jAr{v) = -j[p,jv]rE* +J^[p,v]rE* = Jv, 

Nj{v) = - [p,J'^v]te-' + J[p,Jv\te* = -Jv. 

and using (1.5.23) it results the conclusion. □ 
This connection is induced by J and p. Its local coefficients are given 

by 

•^"^ = \ " ^^^^^^ " P^^'^^^ + C"M^7a) ■ (1-5-26) 

Definition 1.5.7 An adapted tangent structure J on TE* is called homo- 
geneous if 

jO-cJ = -J. 

Notice that J is homogeneous if the local components taisix,^) are 0- 
homogeneous with respect to p. 

Proposition 1.5.8 Let J be a homogeneous adapted tangent structure. A 
section p ofTE* is J— regular if and only if 

Jp = C. 

Proof. If p is J"— regular then 

9 Pa 

is 0-homogeneous, hence must be 1-homogeneous with respect to p, there- 
fore = Pat""^ 1 that is equivalent to Jp = C. Vice versa, if Jp = C then 
= pj'^l^ and thus 

^-t^^ + a^-t^^ 

OjLt-y OPs 

which ends the proof. □ 
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Remark 1.5.6 (i) Based on the above result, the local expression for a J— 
regular section with J a homogeneous adapted tangent structure is 

p = fiat''<^X0+Par''. (1.5.27) 

(ii) The coefficients (1.5.26) generated by p from (1.5.27) can be written 
in the following form 
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1.5.3 Hamilton sections 

Definition 1.5.8 A section on TE* is called a Hamilton section if it is 
J— regular and 

^■4>^E = 0, 

where uje is the canonical symplectic section. 

If tp = ^"Qa + Pa'P'^ then the condition C^oje = is expressed locally 
by [IM] 



^^^~dc^ 1*^/3-^7" ~ ^a^lPj - Pl^ap- 

Now, we deal with some generalizations of the Hamilton sections. 

Definition 1.5.9 The section ip = ^°^Qa + PaP°^ on TE* defines a mechan- 
ical structure if ip is J —regular and 

i^e{Jpi,P2) = uje{JP2,Pi), (1.5.29) 
for any pi,p2 G T{TE*). 

The definition is equivalent with the symmetry of t"'^ = which 

OjJ.ce 

means that the property (1.5.28) i) is fulfilled. 

Proposition 1.5.9 The section ip on TE* defines a mechanical structure 
if and only if 

C^UJE {1^1,1^2) = 0, 

whenever vi,V2 £ TiVTE*). 

Proof. Let us consider pi,p2 € T{TE*) and ui = J^pi, 1^2 = 'Jp2- Then 

C^uje{Jpi,JP2) = 'P{^^e{Jpi,JP2)-^^e{{^'4,Jpi,JP2)- 

-0Je{{JPI.C.^JP2) = U)E{{Mpi,Jp2)+i^E{{Jpl,Mp2)- 
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Considering M pi , Mp2 for pi , p2 in the previous relation and using the prop- 
erties (1.5.23) and = Id we obtain 

C^oje{Jpi,JP2) = C^i^E{J.N'pi,jMp2) = u^e{pi,JP2) + u^e{Jpi,P2), 
which ends the proof. □ 

Definition 1.5.10 The section ip on TE* is a semi-Hamilton section if ip 
is J'— regular and 

iviC^UE) = 0, (1.5.30) 

whenever v G T(yTE*). 

In the case of a semi-Hamilton section on TE* only the conditions 
(1.5.28), i) and ii) are satisfied. 

Let us consider a section if) defining a mechanical structure on TE* , then 
we can find the other connection. Indeed, we consider 

g{pi,P2) = -^e{Jpi,P2), 

for pi,p2 S T{TE*). It results that g{pi,p2) = g(p2,pi) and g{v,p) = 
whenever v G T{VTE*). The local coordinate expression of g is given by 

Proposition 1.5.10 If ip defines a mechanical structure on TE* then the 
section N' defined by 

oje{N'pi,P2) = (£v;ff)(pi,P2), Pi,P2 G TiTE*), (1.5.31) 

determines a connection on TE* . 

Proof. First, we show that the following relation holds 

{C^g){pi,P2) = -{C^uje){Jpi,P2) +uje{Mpi,P2). 
Indeed, we get 

(^v5)(Pi>P2) = il^{g{pi,p2)) - g{C^pi,p2)- g{pi,C^p2) = 

= ~llj{ujEiJpi,P2)) +(^EiJi^ipPl),p2) + W£;(Jpl,/:^p2)• 
But the relation 

tp{l^E{Jpl,p2)) = {I-"iI,OJe){J Pi, P2) + I-^e{I-"iI,{J Pi), P2) +i^E{Jpl,J^iljp2), 

yields 

{C^g){pi,P2) = -{C^loe){Jpi,P2) + we(-(/:^J)pi,P2)- 
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Also, it results 

{C^uje){Jpi,P2) = uje{{M - M')pi,p2). (1.5.32) 

But, 

'^eWJpi,P2) = {C^g){J Pi, P2) = 
= i^igiJ Pi, P2)) - g{{C^J)pi,P2) - g{Jpi,C^p2) = 
= g{Mp\,p2) = -^e{JM P\,p2) = -i^E{Jpi,P2), 
whence M'J = - J. Next 

u;e{JJ^'pi,P2) = -i^E{M'pi,Jp2) = -{C^g){pi,Jp2) = 

{C^0Je){JPI,JP2) - '^e{-N'PI,JP2) = ^e{JM P\,P2) = ^e{J P\-,P2)-, 

whence JM' = J , which ends the proof. □ 
In local coordinates the connection M' has the coefficients given by 

Using (1.5.32) we obtain 

Remark 1.5.7 In the case that ip is both semi-Hamiltonian and mechan- 
ical sections and moreover Hamiltonian, then the connections N and N' 
coincide. 
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1.5.4 Hamiltonian formalism on Lie algebroids 

Let us consider a differentiable and regular Hamiltonian T-L : E* ^ i.e. 
the matrix 

is nondegenerate. 

Any regular Hamiltonian 71 on E* induces a pseudo-Riemannian metric 
on VTE* (the metric tensor is g°'^{q, /x)) therefore, it induces a unique sym- 
metric adapted almost tangent structure (denoted J^u) such that (1.5.18) 
is verified. Moreover, this is a tangent structure i.e., is integrable. 

A ^7— regular section induced by the regular Hamiltonian T-L is 

f)l-/ 

OPa 

Since uje is a symplectic section on the Lie algebroid {TE*, [•, ■]-j-e*,(^^) and 
dn G r{TE*), we get 

Remark 1.5.8 There exists a unique section py^ € T{'TE*) such that 

ip^ujE = d^T-L, 

and p-}{ is a Hamilton section, i.e the relations (1.5.28) are fulfilled. 
With respect to the local basis {Qa,'P'^} , the local expression of p-^ is 



Thus, the vector field a^{p-}i) on E* is given by 



dn d ( , dn dn\ d 



and consequently, the integral curves of py (i.e. the integral curves of the 
vector field cr^{pn)) satisfy the Hamilton equations 

dq^ ,• dn diia i dn ^-y dn ^ 
The Theorem 1.5.2 yields |103j 

Corollary 1.5.1 The symmetric nonlinear connection M = —Cp^J-y has 
the coefficients given by 

1 d^n 
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/ -1/1 _ 1^'^ dgap &H 

2S i/ie Poisson bracket. 

For the particular case of the cotangent bundle see [HTl [901 HZ! 
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1.5.5 Duality between the Lagrangian and Hamiltonian formal- 
ism 

We have a local diffeomorphism <I> from E* to E, locally given by 

x' = q\ y»=e{q,f^) = ^- (1-5.36) 

and its inverse has the following local coordinates expression 

dC 

q' = x\ = (a{x,y) = —, (1.5.37) 

where 

^x,y) = Cay" -^(x,/u), 

and the components (^aix,y) define an 1-section on E. 

From the condition for to be the inverse of $ we get the following 
formulas 

V/3(Ca) o ^ = 9a^, 

= o $-i)V/3, (1.5.38) 

where is tangent map of <I> (see [64J) and g°^'^ = 5^"/9/U^, g"^ gp-y = S^. 

Theorem 1.5.3 Let p be a J'— regular section on TE* induced by the hy- 
perregular Hamiltonian %, and ^ : E* ^ E the global diffeomorpfism given 
by (1.5.36), then the section <I>*p is a semispray on TE whose induced con- 
nection N is the image by <J> of the connection N defined by p. 

Proof. Let us consider p = ^"Qa + Pa'P'^ and we obtain the semispray 

^ = <^,p = e°(A'„ + {Qa{^^) o ^-'W + Paig""^ o $-i)V;3 = y°Af, + l9"V«, 
where t?" is given by 

o CD = e^Q^(r) + pisv^in = Pin- 

Considering $ the map induced by <I> at the level of tensors and using (1.5.38) 
we obtain 
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which yields 

N = -C^J = -U,p^J = -HCpJ) = ^Af, 

thus A'" is the image of by □ 
This theorem shows that the sphtting VTE* © HTE* defined by M is 
mapped by into the sphtting VTE © HTE defined by A^, and it results: 

Corollary 1.5.2 The following equalities hold 

^Jl = 6a, K^6a = 6*a, 

Map{q,l^) = - {Nl{x,y) + Qa{£?))gp^. 

dC 

Proof. We have on one hand 
and on the other hand 

Therefore, we get A^o-y = — {N^^ + Qo(^^)) g^£. Next, since is a diffeo- 
morphism, we obtain 

and using (1.5.38) it results the relation between the curvature tensors of 
connections on TE and M on TE* . By direct computations the other 
relations are obtained. □ 
Let us consider 5 a semispray on TE and :£'—>■ the diffeomor- 
phism given by (1.5.37). Then we set [109j : 

Theorem 1.5.4 The section p = ^^^S is a semi-Hamiltonian section on 
TE* if and only if the nonlinear connection on TE induced by semispray is 
the canonical nonlinear connection induced by the regular Lagrangian. 

Proof. Considering 5 = y°'Xa + S^Va we have from (1.5.38) 

K's = r + i-eg^sQaie) + S'^ga^) v\ 
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The conditions (1.5.28) i), ii) and (1.5.38) lead to 

{XeiU) - A'„(Ce) + Ve{S')gee. + e^eigae) + S^V^igae) + QeLU)+fLi, = 0, 
and using (1.3.25) it follows 

which concludes that (change a with 9 and totalizing) 

Sigae) - Nig.a - Kgee + / {Ll^gea + Ll^gee) = ^^'U^^ Vg (Va, Ve) = 0, 
and 

MCe) - MCa) = Nlg,^ - Nlg,e + CeLle, 
which is equivalent with 



-i d'^L i dL ^ (1.4.20) /•, , X „ 



From the Theorem 1.4.2 we get that the section p = ^^^S is semi- 
Hamiltonian if and only if the connection induced by semispray is the canon- 
ical nonlinear connection (1.4.18). When p = $^^5 is a Hamiltonian section 
we have the same result, but moreover, the condition (1.5.28) Hi) seems to 
lead to the third Helmholtz condition on a Lie algebroid (1.4.22) and the 
work is in progress. □ 
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1.6 Dynamical covariant derivative and metric non-linear 
connection on TE* 

Definition 1.6.1 A map V : T(rs*\{0}) 'Z{TE*\{{)}) is said to be a 
tensor derivation on 7'-E'*\{0} if the following conditions are satisfied: 
i) V is W-linear 

a) V is type preserving, i.e. VT^(r^*\{0}) C 11{TE*\{0}), for each 
(r, s) E N X N 

Hi) V obeys the Leibnitz rule V(i-*(8>5') = VP^ S + P(S>'VS for any tensors 



We consider the M-linear map Vp : r(r£;*\{0}) r{TE*\{0}) by 

v,x = h[p,hx]rE'+v[p,vX]rE', yx G r(r^*\{o}). (i.e.i) 

where /) is a ^7— regular section and it fohows that 



Any tensor derivation on TE*\{{f\ is completely determined by its actions on 
smooth functions and sections on TE*\{Q} (see |117] generalized Willmore's 
theorem, p. 1217). Therefore, there exists a unique tensor derivation Vp on 



We will call the tensor derivation Vp, the dynamical covariant derivative 
induced by the ^/-regular section p and a nonlinear connection N . 

Proposition 1.6.1 The following equations hold 



P,S on TE*\{0} 

iv) V commutes with any contractions. 



Vp(/X) = p{f)X + /VpX, V/ G C^iTE*\{0}), 



(1.6.2) 



TE*\{0} such that 



\c°°iTE*\{o})— P, Vp |r{r£;*\{o})- Vq. 




-r 




(1.6.3) 



n 



(1.6.4) 



where 



(1.6.5) 



The action of Vp on the Berwald basis has the form 




V,51=\i[p,51]te' 
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For a pseudoRiemannian metric g on TE* the action of Vp is given by 

V.giX, Y) = p{g{X, Y)) - g{V Y) - g{X, VpF), (1.6.6) 
which in local coordinates leads to 

gf ■■= Vp5 [V-,V^) = p{g-^)-g^^ (rW,, - (t^W,, - . 

Definition 1.6.2 A nonlinear connection is called metric or compatible 
with the metric tensor g if Vpg = 0, for all J' -regular sections p, that is 

p{g{X, Y)) = g{VpX, Y) + g{X, VpY), VX, Y G T{VTE*). (1.6.7) 
Theorem 1.6.1 The connection M with the coefficients 

M'aP = ^ap + ^taegi3^gj^, (1.6.8) 

is a metric nonlinear connection. 

Proof. Let us consider the dynamical covariant derivative induced by p 
and M given by 

v.. {v,v^) = .b"^) - - ^) - (t^W, - ^) . 

and using (1.6.8) it follows 

(V^V^) = gf - \g'H'^n,r9ep9j - \rtf'^hrgev9f = 0, 

that is is a metric nonlinear connection. □ 
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1.6.1 Nonlinear connection induced by a J -regular section 

If J is an adapted tangent structure and p is a regular section then [50] 

is a nonlinear connection on TE* with local coefficients given by (1.5.26) 



7 



Definition 1.6.3 The Jacobi endomorphism ip is given by 

= v[p,hX]rE'- 

Locally, from (1.6.4) we obtain that ip = IZa^Q" , where Tla/3 is 
given by (1.6.8). TZai3 are the local coefficients of the Jacobi endomorphism. 

Proposition 1.6.2 The following result holds 

ip = ipQ. + v[vp, \iX]rE* ■ 

Proof. Indeed, ipiX) = v[/9, hXjf^;* = v[hp, hXjf^;* + y[yp,\iX]^e* 
and Q.{p, X) = v[hp, \iX]rE* , that is ipiX) = n{p, X) + v[vp, hX]rE* ■ □ 

Remark 1.6.1 If p is a horizontal section p = hp, then we obtain 

ip = ip^, 

and locally it follows 

P-y ~ C -^ay-i 

which yields 

= Tleape- (1-6.9) 

and the local relation between the Jacobi endomorphism and the curvature 
tensor of the nonlinear connection is obtained. 
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1.6.2 Hamiltonian case 



In what follows, we consider a regular Hamiltonian H : E* ^ R, that is the 
matrix 

is nondegenerate. This regular Hamiltonian determines a symmetric canon- 
ical nonlinear connection given by (1.5.35). We have the following result: 

Theorem 1.6.2 The canonical nonlinear connection induced by a regular 
Hamiltonian % is a metric nonlinear connection. 

Proof. Introducing the coefficients (1.5.35) into the expression of the 
dynamical covariant derivative and using (1.5.33) we obtain 



— O", 



d^n , on dg^e ri dn f \ 



Prom the equalities 



f^,<n%^ = -,'^,,,'-f = (1.6.10) 

g^^g-^^-p^ = -g^f^gJ41 = -^4^^ L], = -Ll (1.6.11) 
^ ^ Qqi y Qqi dqi ' eS Se \ ' 

by direct computation, it follows that Vpg (V'^,V^) = 0, which ends the 
proof. □ 

Theorem 1.6.3 The canonical nonlinear connection induced by a regular 
Hamiltonian is the unique metric and symmetric nonlinear connection. 
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Proof. Let us consider a metric and symmetric nonlinear connection 

lave 



with the coefficients N-ye- Then we have 



and using (1.5.33) we obtain 



But the connection is symmetric (1.5.10) and using (1.6.10), (1.6.11) we get 

TT 

idndgej_ _ idndgej_ dT^dge^ _ 

"'dq^ 5w dq^ 5^ " 

n \r , TT i ^ ^ i ^ ^ T a r a 

and we get the coefficients (1.5.35), which ends the proof. □ 
For the particular case of the cotangent bundle see |100t I105j 
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1.7 Poisson-Lie algebroids 

Poisson manifolds were introduced by A. Lichnerowicz in his famous paper 
[66] and their properties were later investigated by A. Weinstein |124j . The 
Poisson manifolds are the smooth manifolds equipped with a Poisson bracket 
on their ring of functions. We remark that the cotangent bundle of a Poisson 
manifold has a natural structure of Lie algebroid. In this chapter we study 
the Poisson structures on Lie algebroids ^3 11U7| . 

1.7.1 Poisson structures on Lie algebroids 

The Schouten-Nijenhuis bracket is given by |120j 

p 1 

[Xl^...^Xp,Yl^...^Y,] = (-if+i ^^(-iy+^'[x„y,-] aXi a ... a 

i=l 3=1 

Xi A... ^Xp^Yl^ ...a Yj a... a 

Let us consider the bivector (i.e. contravariant, skew-symmetric, 2- 
section) H G T{^^E) given by the expression 

Tl=]^^T'^P{x)so^^sp. (L7.1) 

Definition 1.7.1 The bivector H is a Poisson bivector on E if and only if 
[n, n] = 0, where [■,■] is Schouten-Nijenhuis bracket. 

Proposition 1.7.1 Locally the condition [n,n] = is expressed as 

(a,e,6) 

If n is a Poisson bivector then the pair [E, 11) is called the Poisson-Lie 
algebroid. The Poisson bracket on M is given by 

{/l,/2} = n(d^/i,d^/2), /l,/2 G C°°(M) 

We also have the bundle map vr* : E* ^ E defined by 

TT#p = ipU, pGT{E*). 
Let us consider the bracket 

[p, 9]^ = C^^pO - C^#gp - (i^(n(p, 9)), 
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where C is Lie derivative and p, ^ € T{E*'). With respect to this bracket and 
the usual Lie bracket on vector fields, the map a : E* ^ TM given by 

a = a o TT^, 

is a Lie algebra homomorphism 

The bracket [., .]7r satisfies also the Leibniz rule 

[pje]^ = f[p,e]^ + a{p){f)e, 

and it results that {E* , [., .Jtt,?) is a Lie algebroid [61] . 

Next, we can define the contravariant exterior differential (F : /\^{E*) — )• 
/\'+\E*) by 

k+l 

d''uj{si,...,Sk+i) = ^ (-l)*+^5(si)a;(si,...,Si,...,Sfc+i) + 

1=1 

+ ^ {-iy^^Uj{[SiSj]^,Sl,...,Si,...,Sj,...Sk+l)- 

l<i<j<k+l 

Accordingly, we get the cohomology of Lie algebroid E* with the an- 
chor a and the bracket [., .Jtt which generalize the Poisson cohomology of 
Lichnerowicz for Poisson manifolds |120j . 

In the following we deal with the notion of contravariant connection on 
Lie algebroids, which generalize the similar notion on Poisson manifolds [39], 
[66]. 

Definition 1.7.2 If p,9 £ T{E*) and e T{E) then the linear con- 

travariant connection is an application D : r(£^*) x r(£^) — )• T{E) which 
satisfies the relations 

i) Dp+g^ = Dp^ + De^, 

ii) Dp{^ + -^) = Dp^ + Dp^, 
Hi) Dfp^ = fDp^, 

iv) Dp{f^) = fDp^ + a{p){f)^, f G C-(M). 

The contravariant connection induces a contravariant derivative ■ 
T{E) —7- T{E) such that the following equalities are fulfilled 

Dp{fe) = fDpe + a{p){f)e, /gC"»(m), p,eev{E*) 

In the case where the contravariant connection D is induced by a covari- 
ant connection V on a Lie algebroid E (see [30j) we have Dp = V^#p. 
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Definition 1.7.3 The torsion and curvature of linear contravariant con- 
nection are given by 

T{p,e) = Dpe-Dgp-[p,eu 

R{p, 0)p = DpDgp - DgDpP - £'[p,6l]^M, 

where p,9,p, er{E*). 

The curvature tensor satisfies the equaUties 

R{p,e) = -R{e,p), Rifp,e) = fR{p,e). 

The Bianchi identities have the following form 

J2iDpR{0,fi) + R{T{p,e),p)) = O, 

(Rip, 0)p - TiT{p, 9), p) - DpT{e, p)) = 0. 
In local coordinates we define the Christoffel symbols F^^ by the formula 

Dsc^s^ = rf 

and under a change of coordinates 

x*' = x*'(x*), = l,n on M ,^ ^ 

= y", a, a' = l,m on E, 

corresponding to a new base s"^' = s"', these symbols transform according 
to 

rf = «^:^,rf + A'^Ay.^TT'^^. (1.7.4) 

If we denote T{s'^,s^) = T^'^s'^ and i?(s",s^)sT = Rf^s^ then, under a 
change of coordinates, we obtain that 

-^7' — ^7'-'7 ' -"-5' ~ ^7 ^^'-"■(S • 

Proposition 1.7.2 The local components of torsion and curvature of linear 
contravariant connection are 

rpap _ ya/S _ y/Sa _ a-y t J3 , r a _ t ^" 

^£ -^e " -^7£ -^7£ ^£ Q^i ) 
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Definition 1.7.4 A tensor field T on E is called parallel with respect to D 
if and only if DT = 0. 



Definition 1.7.5 A contravariant connection D is called a Poisson connec- 
tion if the Poisson bivector 11 is parallel with respect to D. 

Remark 1.7.1 If the Poisson connection D is induced by a covariant con- 
nection V (i.e. DH = 0, Dp = V^#p, ir'^Dpcl) = V ^#p'K't^4') then the torsion 
and curvature tensors of the both connections are related by the following 
equalities 

T^{7r*p,Tr*0) = Tr*T^{p,e), 
R^{tt*p, Tr*e)Tr*fi = iT*R^{p, 0)p, Vp, e,pe V{E*). 

Let r be a tensor of type (r, s) with the components T^l"']l and 9 = 6as'^ 
a section of E* . The local coordinates expression of contravariant derivative 
is given by 

DgT = e^;;^rsi, ®---®Si^®s^^®---® s^% 

where 

Qr-£il...ir r s 

rpii...ir /a as i ji---js , \ ^ / -piaarpii...e...ir\ \ ^ / -pearpii.-.ir \ 

^h-3sl ~^ Qy.i Z^y-s ^h-3s ) Z^y- 3b-^h -e...js) ' 

o=l 6=1 



and / denote the contravariat derivative operator. 

Let V 
We have 



Let us consider a contravariant connection D with the coefficients T^^ . 



Proposition 1.7.3 The contravariant connection D with the coefficients 
given by 

T°;^ = rf - ^7r^.7r-/^ (1.7.5) 

is a Poisson connection. 

Proof. Considering / the contravariant derivative operator with respect 
to contravariant connection D, we get 



' 2 ' 2 ' 



□ 
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Proposition 1.7.4 a) The functions 

rf-^^, (1.7.6) 

are the coefficients of a contravariant connection. 

h) The contravariant connection with the coefficients given by (1.7.6) is 
a Poisson connection if and only if 

TT^'^vr^^^^^ = 0. (1.7.7) 

a,e,S 

Proof, a) Using the change of coordinates (1.7.3) and the fact that the 
structure function cr^ change by the rule [69] 



A' Aa' _ „i 

we obtain that the coefficients (1.7.6) satisfy the transformation law (1.7.4). 
b) We obtain that 

and using the equahty (1.7.6), it results 

(",/3,7) 

Prom the condition [11, 11] = 0, locally given by the equation (1.7.2), it results 
that 11^"'/°' = if and only if the required relation is fulfilled. □ 

Remark 1.7.2 Considering 

pa/3 ^ i ^ 

in relation (1.7.5) we obtain a Poisson connection D with the coefficients 

-al3 . 97r"^ 1 



which depends only on the Poisson bivector and structural functions of Lie 
algebroid. 

Proposition 1.7.5 The set of Poisson connections on Lie algebroid are 
given by 

7 — 7 ' ■yv e i 

where 

and D{Tj^) is a Poisson connection with X^^ an arbitrary tensor. 
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Proof. By straightforward computation it results 



7]-'" ' / = TT cr, — — ^ hip. TT ' + i- „ W = 

^/37/« + ivr^^X,^" - ^TT^^XJ" + ^vr^'^Xj" - ^vr'^^X,^" = 0, 

because tt^'''/" = 0, which ends the proof. □ 
The image of the anchor map cr{E) C TM defines an integrable smooth 
distribution on M. Therefore, the manifold M is foliated by the integral 
leaves of (t{E), which are called the leaves of the Lie algebroid. A curve 
u : [to)^i] E \s called admissible if a(u{t)) = c(t), where c(t) = iT{u{t)) is 
the base curve on M. It follows that u{t) is admissible if and only if the base 
curve c{t) lies on a leaf of the Lie algebroid whereas two points can be joint 
by an admissible curve if and only if they are situated on the same leaf. We 
can choose a smooth family t — t- 9{t) G E* of 1-form such that a9(t) = c{t). 
We shall call the pair {u(t), 6{t)) the dual curve. 

Definition 1.7.6 Let (u{t), 0{t)) a dual curve on E. We say that (u{t), 
9{t) ) is a geodesic if 

{De9)^^t) = 0. 

In local coordinates we obtain that a curve 

iu{t),9it)) = {x\t),...,x^it),9iit),...,9m{t)) 

is geodesic if and only if it satisfies the following system of differential equa- 
tions 

^ =ai^Tr-^{x\t),...,x-{t))9a{t) 

(L7.8) 

^ = -T^^{x\t),...,x-it))9i,9^. 
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1.7.2 Compatible Poisson structures 

Let us consider the Poisson bivector on Lie algebroids given by the relation 
(1.7.1). We obtain 

Proposition 1.7.6 The complete lift ofH on TE is given by 

2^^;^ - T^'^L^^j y'^Va A Vp. (1.7.9) 
Proof. Using the properties of the vertical and complete lifts we obtain 

= \ t:"^ si A si + \^-^{sl A s} + si A s^) = \ tt"^ V« A V;3+ 
+i7r"/3 ((A-^ - Li^y^Vs) A + A [X^ - L^f.^y')^^ = 

□ 

Proposition 1.7.7 The complete lift W is a Poisson bivector on TE . 

Proof. Using the relation (1.7.9) by straightforward computation we 
obtain [II'^,!!'^] = which ends the proof. □ 

Proposition 1.7.8 The Poisson structure 11^ has the following property 

n'^ = -/:c^^ 

which means that {TE, li^) is a homogeneous Poisson manifold. 
Proof. A direct computation in local coordinates yields 

Ccli' = jOyev, {ir^^Xa AVp + i^i,^ - 7r^'3L^^)yTVa A V^) 

= CysvA^^^^a) A - TT^^Xa AVp + Cyey^lai^^yWa) A 

-H^y^'^c. AVfs- {Cy^vy^L^^yWc,) A + TT^^L^^y^Va A 
= -1^'^^X^ A - ^4^y^Va A + TT^^LI^y-^Va A 

□ 
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Definition 1.7.7 Let us consider a Poisson bivector on E given by (1.7.1) 
then the horizontal lift o/II to TE is the bivector defined by 

Proposition 1.7.9 The horizontal lift is a Poisson bivector if and only 
if n is a Poisson bivector on E and 

n(^iTT*pf,{n*ef'^=o, Mp,eeT{E). (1.7.10) 

Proof. The Poisson condition [11,11] = leads to the relation (1.7.2) 
and [n^,n^] = yields 



(e,(5,a) 



and tt'^^tt'^^TZ'^p^ = 0, which is the local expression of the relation (1.7.10). □ 
We recall that two Poisson structures are compatible if the bivectors loi 
and UJ2 satisfy the condition 

[uJi,U2] = 

By straightforward computation in local coordinates we get: 

Proposition 1.7.10 The Poisson bivector is compatible with the com- 
plete lift if and only if the following relations hold 

vr '^vr (— — ) — vr 'vr L':^. = U, 

^ dy^ dy'' ' '"^ 

where we have denoted 



' dx'' 

See [HI] for the particular case of tangent bundle. 
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1.7.3 Canonical Poisson structure 

On Lie algcbroid (TE*, [•, •vtJt-^'* , cj-'^) we have the canonical symplcctic sec- 
tion ue given by (1.5.6) which induces a vector bundle isomorphism 

i^^^:E*^E, i^cjE e E* ^ C & E. 
Definition 1.7.8 The canonical Poisson bivector is given by 

It follows that 

A{dF, dG) = -ueWF), ^(dG)), i^, G G C~(£;*) 
and in local coordinates we get 



Remark 1.7.3 The Schouten-Nijenhuis bracket [A, A] leads, locally, to the 
expression 

and [A, A] = follows from the structure equations on Lie algebroids (1.2.9). 
Definition 1.7.9 Let us consider a Poisson bivector on E given by 

Ii=^'K^P{x)SahSp 

then the horizontal lift of 11 to TE* is the bivector defined by 

Proposition 1.7.11 The horizontal lift 11^ is a Poisson bivector if and 
only if n is a Poisson bivector on E and 

n(^{7T*pf,{7r*ef'^=o, \/p,eer{E*) (1.7.11) 

Proof. The Poisson condition [W, W] = leads to the relation 

(tt Ptt^ L^^ + tt P(j^^) =0 

(a,£,<5) 
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and [n^, n^] = yields 



(£,6,a) 



and it results 7r"^7r''''^7?./j-y£ = 0, which is the local expression of the condition 
(1.7.11). □ 



Proposition 1.7.12 If the connection J\f on TE* is defined by a linear con- 
nection V with the coefficients T^^^ on the Lie algebroid E, the the bivector 
has the following form 

= l7r"'3Q„AQ^+^7r°^/x^Mer2:,rj57'^AP'^+7r"^Ai7r^.2aA7'^ (1.7.12) 

Proof. The coefficients of the nonlinear connection have the form A^^^ = 
/^■y^aj3 ^^'^ introducing the relation 5^ = Qa + J^apT^^ into the expression 
of it results (1.7.12). 

Proposition 1.7.13 If M is a symmetric nonlinear connection then the 
canonical Poisson bivector has the form 

A = r'^A 61. 

Proof. We have 

= A 5* - iNapV" AV^ - A 

= A 5* - WapV^ KV^ - WapV^ A 

Proposition 1.7.14 // is a Poisson bivector and N is a symmetric 
nonlinear connection, then 11^ is compatible with the canonical Poisson 
structure A if and only if the following relations fulfilled 

7r"^7^a7e = 0. (1.7.14) 
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Proof. If M is symmetric then A/'q./j — AZ/ja = IJ^'yL'l^ji and with respect 
with the basis {^^i^"} it results A = A (5*. By a straightforward com- 
putation we obtain 



+ 



2 V dp^o 

and [n-^, A] = is equivalent with the relations (1.7.13), (1.7.14). □ 



Remark 1.7.4 If the nonlinear connection M is defined by a linear connec- 
tion V with the coefficients F^^ on the Lie algebroid E then we obtain the 
conditions 



,, i ^« "7 I -pe -pe re T^' r f re 1 _ n 
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2 Optimal Control 



The Lie geometric methods in control theory have been appUed by many 
authors (see for instance. [TSl ESI [13 EI])- One of the most important is- 
sues in the geometric approach is the analysis of the solution to the optimal 
control problem as provided by Pontryagin's Maximum Principle; that is, 
the curve c{t) = {x{t),u{t)) is an optimal trajectory if there exists a lift- 
ing of x{t) to the dual space {x{t),p{t)) satisfying the Hamilton equations, 
together with a maximization condition for the Hamiltonian with respect 
to the control variables u{t). The purpose of this part is to study the 
drift less control affine systems (distributional systems) with positive ho- 
mogeneous cost, using the Pontryagin Maximum Principle at the level of a 
Lie algebroid in the case of constant rank of distribution. Author's papers 
Ha ESI EH [Ml [Ml [Ml [M] are used in writting this part. 

We prove that the framework of Lie algebroids is better than cotangent 
bundle in order to solve some problems of drift less control affine systems. 
In the first chapter the known results on the optimal control systems are 
recalled by geometric viewpoint. In the next chapter the distributional sys- 
tems are presented and the relation between the Hamiltonians on E* and 
T*M is given. We investigate the cases of holonomic and nonholonomic 
distributions with constant rank. 

In the holonomic case, we will consider the Lie algebroid being just the 
distribution whereas in the nonholonomic case (i.e., strong bracket gener- 
ating distribution) the Lie algebroid is the tangent bundle with the basis 
given by vectors of distribution completed by the first Lie brackets. In the 
both cases illustrative examples are presented. Also, the case of distribu- 
tion D with non-constant rank is studied and some interesting examples 
are given. In the last chapter we present the intrinsic relation between the 
distributional systems and sub-Riemannian geometry. Thus, the optimal 
trajectory of our distributional systems are the geodesies in the framework 
of sub-Riemannian geometry. We investigate two classical cases: Grusin 
plan and Heisenberg group [39] , equipped with positive homogeneous costs 
(Randers metric 0). 
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2.1 Geometric viewpoint of the optimal control 

Let M be a smooth n-dimensional manifold. We consider the control system 

where x £ M and the control u takes values in an open subset Q of R"^. 
Let xq and xi be two points of M. An optimal control problem consists of 
finding the trajectories of our control system which connects xq and xi and 
minimizing the cost 



mm / L{x{t),u{t))dt, x{0) = xq, x{T) = xi, 
Jo 



where L is the Lagrangian or running cost. 

Necessary conditions for a trajectory to be an extreme arc given by Pon- 
tryagin Maximum Principle. The Hamiltonian reads as 

H{x,p, u) = {p, f{x, u))) - L{x, u), p G T*M, 

while the maximization condition with respect to the control variables u, 
namely 

H{x{t),p{t),u{t)) = max H{x{t),p{t),v), 

V 

leads to 

— - 

du 

The extreme trajectories satisfy the Hamilton equations 

. dH . dH 



From a geometric viewpoint the pair (a;*,M°), where i = l,n and a = l,m, 
can be understood as a local coordinate pair of a manifold E, that is fibered 
over M by the projection tt : E ^ M. The functions f^{x,u) are the 
components of a vector field X = P{x,u)-^ along tt, that is, of a fibered 
mapping X : E ^ TM from the bundle (i?, vr, M) to the tangent bundle 
(TM, r, M) such that toX = vr. The admissible curves of the control system 
are the curves j : I C R ^ E such that 

xm) = |(^(7W)- 
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The optimal control problem consists of obtaining the admissible curves that 
minimize the cost 



mm 



and satisfy certain boundary conditions not to be considered here. 

The Hamiltonian if is a real-valued function defined on the fibered product 

T*M Xm E that is given by 

H{^l,v) = {^i,x{v))-L{v), 

for any (//, v) eT*M Xm E. 

The critical equations follow from asking a vector field Xh defined along a 
map pri : T* M x m E ^ T* M to satisfy the symplectic equations 

ixn^ = dH, 

where uj = dx^ A dpi. Since 

dH = —dx' + —dpi + TT^du'', 
ox^ opi ou"- 

we obtain that the solution of the above equations is the vector field 

dH d dH d 



Xh 



dpi dx^ dx^ dpi ' 
defined on the subset 

l^ = » 

du"- 

of T*M X M E, and therefore, the critical trajectories are the integral curves 
of the above vector field, namely 

dH ^ ■. dH . dH 

= 0, X* = — , pi = - — . (2.1.2) 

au" dpi dx^ 

By a control system on the Lie algebroid (see [71]) ir : E ^ M with the 
control space r : ^ — t- M we mean a section p oi E along r. A trajectory 
of the system p is an integral curve of the vector field (j{p). Given the cost 
function C E C°°(A), we have to minimize the integral of C over the set 
of those system trajectories which satisfy certain boundary conditions. The 
Hamiltonian function % G C°° {E* x m A) is defined by 

'H(/i,'u) = {p,p{u)) - C{u), 

whereas the associated Hamiltonian control system pH is given by the sym- 
plectic equation 
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In local coordinates, the solution of the previous equation reads as 
on the subset where 

^ = 0. 

Therefore, the critical trajectories are given by 

d^-^' 'dt - "'^dji:' 'df-'^^'^W^'' ''^diT^' ^ ^ ^ 
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2.2 Distributional systems 

Let M be a smooth n-dimensional manifold. We consider the distributional 
system (drift less control-affine system) 

m 

x = Y,UiXi{x), (2.2.1) 

i=l 

where x E M , Xi,X2, ...,Xm are smooth vector fields on M and the control 
u = {ui,U2, ■■■,Um) takes values in an open subset 0, of i?™. The vector 
fields Xi,i = l,m, generate a distribution D C TM such that the rank 
of D is constant. Let xq and xi be two points of M. An optimal control 
problem consists of finding those trajectories of the distributional system 
which connect xq and xi, while minimizing the cost 

min / F{u{t))dt, (2.2.2) 
«(•) J I 

where is a Minkowski norm (positive homogeneous) on D. 

Remark 2.2.1 We can associate to any positive homogeneous cost T on 
the Lie algebroid E a cost F on Ima C TM defined by 

F{v) = {T{u)\ue E^, a{u) = v}, (2.2.3) 

where v G {Ima)x C T^M, x G M. 

A piecewise smooth curve c : / C i? — t- M is called horizontal if the 
tangent vectors are in D, i.e. c{t) G -Dc(t) C TM for almost every i G /. Let 
u : I ^ E he an admissible curve projected by tt onto the horizontal curve 
c: I M. The length of the horizontal curve c is defined by 

length{c)= j F{u{t))dt = J F{c{t))dt, 

and the distance is given by d{a,b) = inf length{c) where the infimum is 

taken over all the horizontal curves connecting a and b. The distance is 
infinite if there is no admissible curve that connects these two points. 

Remark 2.2.2 The energy of a horizontal curve is 

E{c) = \jF\c{t))dt, 

and it can easily be proved that if a curve is parameterized to a constant 
speed, then it minimizes the length integral if and only if it minimizes the 
energy integral. 
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For the 2-homogeneous Lagrangian L = and L = we have 

C — L o a. 

Further, L on TM is a Lagrangian with constraints. According to the Pon- 
tryagin Maximum Principle, the Hamiltonian is recast as 

H{x,p,u) = {p, x) — L{x,u). 

If the equations 

dH{x,p, u) 



du 



0, 



permit us to find in a unique way u as a, smooth function of then we 

can write the Hamiltonian system without any dependence on the control. 
This nice situation happens always for distributional systems with quadratic 
cost 



- i=l 



If the cost is not quadratic, then we cannot guarantee that the Hamiltonian 
can be calculated without dependence on the control. However, there exist 
several situations when the Hamiltonian can still be found. 

Proposition 2.2.1 The relation between the Hamiltonian H on cotangent 
bundle T*M and the Hamiltonian % on dual bundle E* is given by 

H{p) = n{fi), fi = a*{p), peT^M, ^eE:. (2.2.4) 

Proof. The Fenchel-Legendre dual of Lagrangian L is the Hamiltonian 
H given by 

H{p) = sup {{p, v) - L{v)} = swp{{p,v) - £.{uy,a{u) =v} 

V V 

= sup {{p,a{u)) - jCiu)} = snp{{a*{p),u)-C{u)} = n{a\p)), 

u 



and we get 

H{p) = n{^i),^i = a\p),p G T;M,/x G El 

or locally 

Ma = (^Ipu (2.2.5) 
where the Hamiltonian Hij)) is degenerate on Kera* C T*M. □ 
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2.2.1 Holonomic distribution 



We assume for the beginning that the distribution D = {Xi, X2, X^) 
is holonomic, which means that S D for every i,j = l,m, i 7^ j. 

In order to apply the theory of Lie algebroids we consider E = D with 
the inclusion as anchor a : E ^ TM. From the Frobenius theorem, the 
distribution D is integrable, it determines a foliation on M and two points 
can be joined if and only if they are situated on the same leaf. 
We consider the following distributional system with positive homogeneous 
cost [SS]: 

X = U^Xi + u'^X2, X = { x'^ 



X 





















^0^ 







and 



rT 

min / F{u{t))dt, F{u) = ^/ {u^f + {v?Y + < e < 1. 

«{■) Jo 

We are looking for the trajectories starting from the point (1,1,0)* and 
parameterized by arclength. The associated distribution D = {Xi,X2) is 
holonomic, because 

^ 9 2 ^ 9 
^1 = ^' ^2 = x — ^ + x — ^ + — ^, 
ax^ ax^ ax^ ax-' 

and therefore [Xi,X2] = Xi. In the case of the Lie algebroid, we consider 
E = {Xi,X2) and the anchor a : E ^ TM'^ has the components 



(2.2.6) 




and we get the Lagrangian 

Using |l8],[83] we can find the Hamiltonian on E* given by 

1 - £2 




(2.2.7) 
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Remark 2.2.3 Using (2.2.4) we can calculate the Hamiltonian H on T*M 
given by H{x,p) = /x = C7*(j)), where 



/Lil \ _ / 1 



Pi 

P2 
P3 



We get that 

TT( ^ ^ I I jPlY , {PlX^ + P2X^ + P^f epi \ rnni>\ 

= 2 l^v 0^ ^ j ■ ^ ^ ^ ^ 

Unfortunately, with H{x,p) from (2.2.8) the Hamilton equations on T*M 
lead to a complicated system of implicit differential equations. 

We will use the geometric model of a Lie algebroid. From the relation 
= L^^i^X^ we obtain the non-zero components = 1; -^21 — ~^ 
while from (2.1.3) we deduce that 

x^ = — + x^—,x^ = x^ — ,x^ = —, 2.2.9 

on . on 

where 



a/Xi (l-£2)2 l_e2 n -p213./JmI)LT^ 

(2.2.10) 



(i - £ ) Y (11:^2)7 + 



The form of the last relations leads to the following change of variables 
Hi(t) = {l-s'^)r{t)seche{t), 



Ii2{t) = a/1 - £^r{t) tanh6'(i). 
In these circumstances we have 



~ ^ _9 V \i 



Y (I - £2)2 ' l-e 

whereas 
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yields 

^J\-£^ -^tanh6'^ = r(-tanh6' + esec/i6'tanh6'). (2.2.11) 

Prom 

we get 

Vl-e^ (^^ tanh^ + 9 sech^^^ = r((l + £)2sech26' - esech^ - esech^^). 

(2.2.12) 

Now, reducing 9 and ^ from the relations (2.2.11) and (2.2.12), we obtain 
\pV--^f = r^£secli^tanli0(£secli0 — 1), 

and 

Vl - e'^9 = r{esech9 - if. (2.2.13) 
The last two relations lead to 

r resech^ tanh 9 
9 £sech0 — 1 ' 

and respectively to 

1 ^ esech^ tanh 9 
r £sech0 — 1 ' 

with the solution 

In \r\ = — liL{e sec h6 — 1) — Inc. 

Therefore 

_ 1 



c (esech^ — 1) 

Since the geodesies are parameterized by arclength, the conclusion corre- 
sponds exactly to the 1/2 level of the Hamiltonian and so we have 

n = ^{l-esech9f = ^. 

Now, c = ±1 and 

1 



r = ±- 



' esech.9 — 1 
From the relation (2.2.13) we have 



dt 1 — esech^ ' 
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and respectively 
The relation 



implies that 



, ci(l — esech6') „ 



Since we are looking for the trajectories starting from the point (1, 1,0)* , 
we have 



and so 



The relation 



leads to 



In = ^ Cl = 1 + £, 
1 + £ 



o,„, , 1 — esech^ , cosh^ — e 
x'^{e) = In r — — = In 



(1 — £)sech0 1 — £ 



^ = -in 

1 ' 



_ C2(l -esech6') 



(1 - £2)seche ' 

whereas from x^{^) = 1 we get C2 = 1 + e. These lead to 

9,^, cosh 6* — £ 

x\e) = 

We obtain also that 

H2 = ^il\x -X I = IJ.IX + x fll, 

and, consequently, 112 = fJ'ix^ + C3. Further, 

^ sinh^ 03(1 — £sech^) 

" (l-£2) sech^' 

Prom x^ (0) = 1 we obtain that C3 = 1 + £ and this yields 
1,^, sinh^ cosh0 — £ 

Remark 2.2.4 If e = we regain the case of distributional systems with 
quadratic cost with the solution 

x^{t) = sinhf + cosht, x'^{t) = cosht, x^{t) = In cosh 
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2.2.2 Nonholonomic distribution 



We assume that the distribution D = {Xi,X2, ...,Xm) is nonholonomic and 
is strong bracket generating [118j . i.e. sections of D and first iterated brack- 
ets span the entire tangent space TM. By a well-known theorem of Chow, 
the system is controllable, that is any two points are connected through a 
horizontal curve (M is assumed to be connected). We also suppose that the 
vectors B = {Xi,X2, ...,Xm, [Xi,Xj]} determine a base in TM. The space 
E = TM with the base is a Lie algebroid over M with at least one struc- 
tural function nonzero. The anchor a : E ^ TAI is just the identity and 
the matrix corresponding to a is determined by the base vectors. 
The control system can be written as 

m 

X = ^ UiXi{x) + OXm+1 + ■■■ OXn, 
i=l 

with 



min / C{u{t))dt. 
«(■) JO 



To solve this minimization problem we consider the Lagrangian 

n 

C = C+ Afc^t^ 

fc=m+l 



(Afc are the Lagrange multipliers) and still work via the maximum principle 

dC 



but at the level of Lie algebroids. We set /ij = 4^ and 



T-L = ^ /ijU* — C, 



i=l 

thus obtaining fij = /xfc = Xk with j = l,m and k = m + l,n. Since C 
is 2- homogeneous with respect to u*, i = 1, m, we get 



i=l 



with the constrains /ifc = A^, A; = m -|- 1, n. 

We consider the following distributional system with positive homoge- 
neous cost [98l: 



X 



U^Xi+V?X2, 



X 





0^ 
















^ j 
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and 



in / J^(u{t))dt, J^{u) = ^J{u^Y + {v?Y + eu^, < £ < 1. 
. ) Jo 



mm 

We are looking for the trajectories starting from the origin and parameter- 
ized by arclength. We have 

V 9 ^ d ^ d 
ox'- OX'^ OX"^ 

d 

X3 = [Xi,X2] = 

and (Xi, X2, X3) = TR^, hence the distribution D =< Xi,X2 > of constant 
rank is strong bracket generating. 

Remark 2.2.5 We can work, as in classical case, directly on the cotangent 
bundle by computing the Hamiltonian H(x,p) = T-L{x,fi) , fj, = cr*{p). Since 



^1 

M2 




we obtain 




epi 



l-e2 

Unfortunately, from the Hamilton equations on T*M a very complicated 
system of implicit differential equations is obtained. 

We will use a different approach. Let us take M = and E = TM 
with the basis {Xi^X2,Xs}. E is a Lie algebroid over M with at least one 
structural function nonzero. The anchor a : E ^ TM is just the identity 
and the matrix of a with respect to the basis of E and TM basis is 

1 
1 
1 



Prom the structure equations of Lie algebroids and the relation [Xq,,X^] 
-7 



L'^aX^ we obtain the non-zero structural functions = 1) -^21 — Now, 
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using the Hamilton equations on Lie algebroids (2.1.3) we get the following 
systems of differential equations 



dn .2 on 

X = , X — 



X = X 



dn2 



(2.2.14) 



and 



r 3 01-L \ dl-L 
fis = ^ \ = ct, 



(2.2.15) 



where 



dn 

dm 

dn 

dn2 



(l-£2)2 



l-e2 



1-e^ 



(l-e2)V(fe^ + & 



We may use to the following transformations 

H^{t) = (1 - e'^)r{t){acosAe{t) - bsmAe{t)), 

fi2{t) = Vl - e'^r{t){asmAe{t) + 6 cos Ae{t)), 
with + = 1. We also have 



yields 



+ 



(M2)2 



(l-£^)2 ^ (l-£2) 



(2.2.16) 

(2.2.16') 

r|. Further, (2.2.15) 



(2.2.17) 



and 



ci (a sin A9 + b cos ^6*) + A9{a cos ^6* - 6 sin 
(1 + £2) (a cos ^0 - b sin A9) - e(l + {a cos A9 - b sin A9f) 

ci ( ^ (a cos A9-b sin A^) - A9{a sin + 6 cos 



= - (a sin ^6* + 6 cos 7461) + £(a cos A9 - b sin ^6*) (o sin A9 + b cos ^6*) , 

(2.2.18) 

where a = Oz^h/El!. Reducing 9 and ^ from (2.2.17) and (2.2.18), we 
get 



ci- =£ (a sin ^6' + 6 cos ^6') (eo cos ^6*- £6 sin ^6'- 1) (2.2.19) 
r 
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and 

ciAe = {l- £{a cos Ae-b sin Ae)f. (2.2.20) 

These lead to 

{l-e^)VT^A f de 



A J {l-e{acosAe-hsuiAe)f' 

Prom (2.2.19) and (2.2.20) we obtain 

dr _ Ae{a sin A6 + h COS. A9) 
r £(acosyl^ — 5sin^0) — 1 

and 

1 



d9 



r = 



c{l - £ {acos A9 - bsinAe))' 

Since the geodesies are parameterized by arclength this corresponds exactly 
to the 1/2 level of the Hamiltonian and we have 

o 1 

'H= — {1-e {acos Ae- bsinAe) f = 
So, c = ±1 and 

1 

r = ± 



1 - e{acosA9 - 6sin^6')' 
Prom (2.2.16) we obtain 

_ (l-£2)(acos^6l-6sin^6') 
~ l-e{acosA9-bsmA9) 



_ Vl - £^ (a sin A9-i-b cos A9) 
~ l-e{acosA9-bsmA9) ' 

Since /12 = Ai;^, we also have x^{9) = ^ — ai. As we are looking for geodesies 
with start from the origin, we have ai = x{i^l^a) therefore 



1 y/1 - e"^ {a sin A9 + b cos A9) ^l-eH 
^ ^ ^~ X {1 - e {a cos A9-b sin A9)) ~ A(l -ea)' 

Prom fii = —\x^ wc get 

2 (1- e'^){b sin A9 - a cos A9) {1 ^ e^)a 
^ ^ ' ~ X{l-e{acosA9 - bsinAO)) ^ X{1 - ea)' 

Pinally, because i'^ = .x^.t^ a straightforward computation leads to 



{l-e'^)Vl^A f {a sin A9 + b cos A9f 
[l-e'^)Ab f a sin 746' + 6 cos ^6' 



A2 y (l-e(acosA^-6sin^e))^ 

d9 



>? J {l-e{acosAe -bsinAe))' 



99 



Remark 2.2.6 Fore = we obtain the distributional systems with quadratic 
cost with the solution ■ 

1,, a sin \t — b(l — cos X9) o ,9 

x^{t) = i ^, a^ + b^ = l. 

A 

2,. b sin Xt + a{l — COS Xt) 

^ V') ~ ) 

o , , t - . „ , ab 2 ^ O'b . b"^ . ^ 

X (t) = — H — sm 2Xt — -T7 cos At + t-^t cos Xt ^ sm Xt. 

^ ' 2A 4A2 A2 y? 
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2.2.3 Distributional systems with no constant rank of distribu- 
tion 

Let us consider in the three dimensional space ii^ the drift less control affine 
system p,02] 

x{t) = u^Xi + u'^X2 + u^Xs, (2.2.21) 

with 

and minimizing the cost 

min j F{u{t))dt, (2.2.22) 

















X 










K 1 







where F = y (u^)'^ + {u'^)'^ + (u^)^ + £u^, < e < 1 is the positive homoge- 
neous cost (Randers metric) . The distribution D is generated by the vectors 
Xi, X2, X3 and we can write D = {Xi, X2, X3}. We observe that 



rankD 



3 if X / 
1 if x = 



In the canonical base (^4,^) of we have X, = f , X2 = x^, 
~ '^'Bz ^^'^ brackets are given by 

[Xi,X2\ = ^ = X^iD, [Xi, X3] = ^ = X5 ^ D, [X2, X3] = 0. 
oy oz 

It results that the distribution is nonholonomic, but is bracket generating, 
because the vector fields {Xi,X2,X^,Xi = [Xi,X2\ ,^5 = [XijXs]} gener- 
ate the entire space E?. 
Prom (2.2.21) we obtain 

dx I I dy 2 2 '^^ 3 3 
— = u = s , — = u X = s , — = u X = s . 
dt dt dt 

The cost function can be written in the form (x 7^ 0) 



2^2 (^3)2 



x"^ X^ 



i=l 
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(Einstein's summation) where = e, 6^ = 0, 6^ = and 

/ 1 

gij = 

V 

The Lagrangian has the form L = and using [83] (Th. 4.5 pp. 191) we 



2 

obtain the Hamiltonian in the form 



1 

H = 

where 



H = -[\/g^jPiPJ-b'p^) , (2.2.23) 



and (7*-^ is the inverse of the matrix gij. In these conditions we obtain that 

1 - 

1 

g^O = \ 





and and it results 



V 



From (2.2.23) we obtain 



or, in the equivalent form 



2(1-^2)2 + 2(1 -£2) l_£2y (l_£2)2 + l _ ^2 ' 

In the case x = we obtain 

^ ^ 1^2 ^ {1 + eful 
2 2 

with the constraints y= 0, 2;= 0. Using Lagrange multipliers we obtain 

Li= L + XiV +A2 z, 
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and from Legendre transformation it results 



,2 



2(l + £)2- 

For x = from (2.2.24) we have 



H 



I ( Pi ^Pl \ _^ Pi 



2\l-e^ 1-ey 2(i + £)2' 
which leads to the equahty 



H\^=o = Hi. 



Next, if we denote 



p,_ pI {pI+pD^^ 

(l-£2)2+ i_e2 ' 

then the Hamilton's equations (2.1.1) lead to the following differential equa- 
tions 

dx _ dH _ {1 + e^)pi £ ./z; 1 (0 0^^ 

dt- dp,- l-£2^^ (l-£2)3Ve' ^ ' 



dy dH P2X epip2X 1 



dt dp2 l-e2 (l-£2)2^' 

dz dH Psx"^ epipsx'^ 1 

dt^d^^ 1 -£2 " (1 -£2)2 



(2.2.26) 
(2.2.27) 



dpi _ dH _ (pI+pI) X ^ epi jpl+pf) X 1 

dt - dx- l-£2 + (l_,2)2 

P2 = o = const. 
D P3 = 6 = const. 

p2 /•„2_|_ 2^2 

In these conditions the relation 6 = ^i_^'2y2 + \_g2 — leads to the following 
change of variables: 



dx 


l-£2 


dp2 


dH 


ill ~ 


dy 


dp3 


dH ^ 


ill ~ 


dz 



x{t) = , . ,„ Pi{t) = (1 - £^) rit) COS Ae{t). 

Va + 

It results G = r'^{t) and from (2.2.25) we get 

dx (l + £2)rcosA0 er ercos^AB 

dt ^ 1 -£2 l-£2 l-£2 ■ 
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But 

dx V 1 — 6^ / - • \ 

— = ^^=^ (rsinAe + rA cosAe] , 
at Va^ + 62 V y 

and it results 

ci (r sin AO + rA0 cos Ae^ = {1 + e'^ycosAO - £r{l + cos^ AO), (2.2.29) 
where we have denoted 

The equation (2.2.28) yields 

dpi y/aPTW . An , An ■ An\ 

—-— = — (— r sm A0 + er cos sm ^t') . 

at -y/i _ £2 

But 

= (l - £^) cos ^6* - ryl ^ sin ^461^ 



dpi 
dt 

which leads to 



ci (r cos AO -rAe sin ^6*^ = -r sin A0 + er cos ^16' sin A0. (2.2.30) 

The equation (2.2.29) multiplied by cos^0, minus equation (2.2.30) multi- 
plied by sin A0 lead to the equation 

ciA^ = {I -e COS A0f , (2.2.31) 

and it results 



I 

J (1 — £CO 



Va2 + 62 j (l-£cosA6') 

Moreover, the equation (2.2.29) multiplied by sin^0, plus equation (2.2.30) 
multiplied by cos A0 lead to the equation 

dv 

ci— = ersi-aAe{ecosAe-l). (2.2.32) 

Prom the equations (2.2.31) and (2.2.32) it results 

dr , zii'mAO 
— = -A- 



r 1 — £ cos AO 

which leads to the following result 

Inr = — lnc(l — ecos^d^?), c e R, 
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and we get 



rm = — -— . (2.2.33) 

c{l -e cos Ae(t)) ^ ' 



Using (2.2.33) the Hamiltonian become 

H = ^(l + £^)r^cos^^6' + ^r^sin^A6'-£r^cosA6' 

= -r^ + \e^r'^ cos^ Ad — er^ cos AQ 
2 2 

= y(l-£COS^6') 

1 

2^ 

Considering the integral curves parameterized by arclength, that corre- 
sponds to fix the level \ of the Hamiltonian, we have c = ±1 and it results 

r = ± ^- — . (2.2.34) 

l-£cosA6' ^ ' 

In these conditions, from (2.2.34) we obtain 



^W = ±^S-^^^i^. (2.2.35) 
V^a2 + i,2 i-£cos^0(i) ^ ^ 

The differential equation (2.2.26) yields 

dy ar^sm^AQ ,^ a sin^ ^6* „ „„x 

-1- = ^ -3— (l-£COsA6') = ^5 -r — . (2.2.36) 

Prom (2.2.31) and (2.2.36) it results 

aA(l-£^)^/' r sin^ AO 
(a2 + 62)^/2 J {I -e cos AO) 

In the same way we obtain 



io? + h^Y' J (1 - £ cos A9Y 



(a2 + 62)3/2 j (i_ ecos AOf ^' ^ ' 

In the particular case of £ = we obtain a distributional system with 
quadratic cost 

F = sj{u^f + {V?f + (^^3)2^ 

with the solution 

, , , sinai 

x{t) = ± , 

a 
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at asin2at 
" 2^ " 4a3/2 ' 

bt asm2at 
^2^~ 4a3/2 ■ 

where a = y/a? + h'^, which are the geodesies in the framework of the so 
cahed sub-Riemannian geometry [14] . 
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2.2.4 Distributional systems with degenerate cost 

Let us consider in the two dimensional space the drift less control affine 
system jl06] 

x(t) = u^Xi + u^X2, (2.2.39) 



with 



and minimizing the cost 



min / F{u{t))dt, (2.2.40) 
"(•) Ji 



where F = \J (u-*-)^ + (n^)^ + is the positive homogeneous cost (Kropina 
metric). 

I have to remark that in the case u < 0, = it resuhs F = that 
is we obtain a degenerate cost (metric). The distribution D is generated by 
the vectors X\ , X2 and we have 



rankD 



2 if x/0, 
1 if X = 0. 



In the canonical base of given by | ^ ) ^ | we can write 



^1 = ^2 = X — 

OX oy 
and the Lie brackets are given by 

[X,,X2] = ^=X,iD. 

It results that the distribution is nonholonomic, but is bracket generating, 
because the vectors fields {Xi,X2,X3 = [Xi,X2]} generate the entire space 
Prom (2.2.39) we obtain 

dx 1 1 

— = u = s . 
dt 

dy 2 2 

— = U X = s . 

dt 

The cost function can be written in the following form (x 7^ 0) 

F= ^(ni)2 + (n2)2 + ui = ^(51)2 + ^ + 5^ = ^gijs'si +Y.b's' 
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(Einstein's summation, i,j = 1,2) where 6^ = 1, 6^ = and 

- M 

- V 1 l/x^ 

The Lagrangian function has the form L = ^F^ and using [83j we obtain 
the Hamiltonian in the form 

H='-('-^)\ (2.2.41) 

where 

is the inverse of the matrix gij. In these conditions we obtain that 



2 V 2pi 
or, in the equivalent form 

2 2 2 

Pi P2 4 P2 2 

8 8pf 4 

The Hamilton's equations (2.1.1) lead to the following differential equations 

dx dH pi P2X^ 



dt dpi 4 4pl 



(2.2.42) 



dy_dH_^^p^ (2.2.43) 



and it results 



dt dp2 2pf 

dpi dH P2X^ p\x 

~dt ~ ~~dx ~ 2pf 2~ 

dp2 dH 

—— = — = U p2 = a = const. 

dt dy 

dx pf — a^x'^ 



(2.2.44) 



dt 4pf ' 

dy a?x^ + apfx"^ 
It " 2^ ' 

dpi a^x a'^x^ 
Ik ~ 2 2^' 
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In these conditions the expression of the Hamiltonian leads to the following 
change of variables: 

^^^^^ r{t)smAeit) ^ =r(i)cosyie(t) (2.2.45) 

We obtain the equations 

dx _ r {cos^ Ae - sin^ Ae) 
'dt ~ 4cos3 A9 ' 

and 

dpi ar sin A6 ar svo? AO 
~dt ~ 2 2cos2 ^6* ■ 

But on the other hand 

dx . rAe .„ 

— = - siuAO -\ cos AO 

dt a a 

=r cos AO -rAe sin AO 

dt 



and it follows 



^ r . ,A cosM6l - sin"^ „ 

- -sinA6' + ^6'cosA6' = , (2.2.46) 

air j 4cos3^6' ^ ' 

1 / A . svnAO siv? AO „ 

- -cos^e-^^sinA^ = — -. 2.2.47 

a\r 12 2cos2^6' ^ ^ 

The equation (2.2.46) multiplied by cos^^, minus equation (2.2.47) multi- 
plied by sin AO leads to the equation 

AdO _ 1 

'^di ~ 4 cos^ AO' ^ ^ ^ 



which yields 



and it follows 



t = — I cos^ AOdO, 
a 



t = - (sin 2A0 + 2A0) 

a 



Moreover, the equation (2.2.46) multiplied by sin^l^, plus equation (2.2.47) 
multiplied by cos AO leads to the equation 

1 r sin / , X 

— = —A TTZ- 2.2.49 

ar 4 cos^ AO 
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Using the equations (2.2.48) and (2.2.49) we obtain 

dr . sin A9 
r cos Ai) 

which leads to the following result 

In |r| + Inci = In |cos A^l , ci G i?"*", 

and we get 

r(t) = — cos^6l(i). (2.2.50) 
ci 

Using the change of variables (2.2.45) the Hamiltonian become 

cos^ AQ sin^ AQ sin^ AQ 
H= + - 7r^ + 



Scos'^Ae 4 8cos2^0' 

and from (2.2.50) we get 

H = ^. 

Considering the integral curves parameterized by arclength, that corre- 
sponds to fix the level 1/2 of the Hamiltonian, we have ci = ±1/2 and 
it results 

r = ±2cosA6'. (2.2.51) 
which together with (2.2.45) lead to the result 

sin 2 Ae(t) 

x{t) = ± ^— ^ (2.2.52) 

The differential equation (2.2.43) yields by direct computation to 

dy _ 2sin2^6' 
dt a ' 

and from (2.2.48) and (2.2.53) it follows 

2A 

d 

which yields 



(2.2.53) 



dy = ^sm^{2Ae)de, 



Finally, the solution is 

.it) = 

„m - -4g(t) smAAejt) 

y^^' ~ a? 4a2 ■ 
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2.3 Sub-Riemannian geometry 

If M is a smooth n-dimensional manifold then a sub-Riemannian structure 
on M is a pair (D,g), where D is a distribution of rank m and 5 is a 
Riemannian metric on D. A sub-Riemannian manifold (M, D, g) is a smooth 
manifold M equipped with a sub-Riemannian structure [1181 185t I14j. 
A piecewise smooth curve c : / C — )• M is called horizontal if its tangent 
vectors are in D, i.e. c{t) G ^c(t) TM, for almost every t G I. In sub- 
Riemannian geometry the length of a horizontal curve c is defined by 

L{c) = l^^/g{c{t))dt, (2.3.1) 

where 5 is a Riemannian metric on D. The distance from a to 6 is 

d{a,b) = mi{L{c)), 

where the infimum is taken over all horizontal curves connecting a to b. 
The distance is assumed to be infinite if there is no horizontal curve that 
connects these two points. 

If locally the distribution D of rank m is generated by Xi, i = l,m, a 
sub-Riemannian structure on M is locally given by a control system 

m 

x = ^Uiit)Xi{x), (2.3.2) 
1=1 

of constant rank m, with the controls u{.). The controlled paths are obtained 
by integrating the above system. If D is assumed to be bracket generating, 
i.e. sections of D and iterated brackets span the entire tangent space TM, 
by a well-known theorem of Chow [23] the system (2.3.2) is controllable, that 
is for any two points a and b there exists a horizontal curve which connects 
these points (M is assumed to be connected). 

The concept of the sub-Riemannian geometry can be extended to a more 
general setting, by replacing the Riemannian metric with a Finslerian one. 
For the theory of optimal control this extension is equivalent to the change 
of the quadratic cost of a control affine system with a positive homogeneous 
cost. Also, the case when the rank of D is not constant may produce inter- 
esting examples. We do not intend to develop a comprehensive study of the 
sub-Riemannian geometry. In the present section we introduce two partic- 
ular sub-Finslerian geometries |l9]: the Grushin plane and the Heisenberg 
group, endowed with some special Randers metrics [3- The geodesies of 
these geometries are obtained by using two different approaches: a direct 
application of the Pontryagin Maximum Principle for the Grushin plane and 
the same principle but combined with some results on Lie algebroids for the 
Heisenberg group. 
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Definition 2.3.1 A sub-Finslerian structure on M is a triple {E,a,J-') 

where 

1) (E,Tr,M) is a vector bundle over M, with the projection map tt. 

2) a : E TM is a morphism of vector bundles. 

3) F is a Finsler metric on E, i.e. T : E [0, oo) and satisfies the follow- 
ing properties: 

a) F is on E\ {0}. 

b) T{\u) = \J^{u) for X> and u e Ex,x E M. 

c) For each y G Ex\{Q} the quadratic form 

9yiu, v) = 2d^^y + + tv}s,t=o 
u,v E Ex, X E M , is positive definite. 

Example 2.3.1 (i) E = M x M™, {Xi, ...,Xm} is a system of m vector 
fields on M, a : E ^ TM, given by 

m 

a{x,u) = ^Ui{t)Xi{x), (2.3.3) 

i=l 

and F is a Minkowski norm on . 

(a) E = D, a : D TM the inclusion and T a Finsler metric on D. 

We can associate to any sub-Finslerian structure {E, a, F) a Finsler met- 
ric on Ima C TM, defined as following 

F{v) = inf{J"(tt)| u e Ex, a{u) = v}. (2.3.4) 

u 

for each 

V e {Ima)x C T^M, x e M. 

Definition 2.3.2 A curve u : I ^ E is called admissible if there is an 
absolute continuous curve c : I ^ M , called horizontal, such that Tr{u{t)) = 
c{t) and cr{{u{t)) = c{t) , t e I. 

The length of an absolutely continuous horizontal curve c{t), t G / is 

length{c) = T{u{t))dt = F{c{t))dt. (2.3.5) 

We can also consider the sub-Finslerian distance 

d{a,b) = ini length{c), 

where the infimum is taken over all horizontal curves connecting a and 6. 
This distance is infinite if there is no admissible curve joining a and b. 
However, if wc assume that Im,a is bracket generating, Chow's theorem 
guarantees that the sub-Finslcrian distance between points is finite. 
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Definition 2.3.3 A length minimizing geodesic or shortly a minimizer is 
an absolutely continuous horizontal curve on M that makes the distance 
between two points. 

Remark 2.3.1 The energy of a horizontal curve is 

E{c) = \jFHc{t))dt, 

and it can easily be proved that if a curve is parameterized to a constant 
speed, then it minimize the length integral if and only if it minimize the 
energy integral. 

If we take the Lagrangians L = and C = ^J^^ we have C = L o a. 
The Fenchel-Legendre dual of L is the Hamiltonian 

Hip) = sup{(p,u) - L{v)] = 



sup < {p, v) + sup {—C{u);a{u) = v} 

V y u 

sup{{p,v) - C{u);a{u) = v} = 

u,v 

sup {{p,a{u)) - £(n)} = 

u 

sup{{a''ip),u)-Ciu)} = nia''ip)) 



Hence 



H{p)=n{fi), fi = cT*{p) (2.3.6) 

p G T*M, n £ E*. The Hamiltonian H on T*M is degenerate on Kera*. 

The Hamiltonian H generates a system of differential equations which 
can be written in terms of canonical coordinates {x,p) in the standard form: 

^, P.= (2.3.7) 

dpi ' * 5x* 

Based on the Pontryagin Maximum Principle we can prove: 

Theorem 2.3.1 Let {x{t),p{t)) be a solution of the Hamiltomr equations. 
Then every sufficient short subarc of x{t) is a minimizing sub-Finslerian 
geodesic. This subarc is the unique minimizer joining its end points. 

Definition 2.3.4 The projected curve x(t) will be called normal geodesic or 
shortly geodesic. 

Remark 2.3.2 Contrary to Finslerian geometry, not every minimizer is 
the projection of a solution of (2.3.7) i.e. is a normal geodesic. Those 
minimizers that are not normal geodesies are called singular geodesies !1118\ 
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2.3.1 Grushin case 

We consider the following distributional system with positive homogeneous 
cost (Grushin Plane) ^49j 



min / T{u{t))dt, T{u 
«(■) Jo 



) = \\u\\+{b,u) , b= {e,Oy, u = {u\u^y, < e < 1. 

x(0) = 0, x(T) = XT- 

We are looking for the geodesies starting from the origin and parameterized 
by arclength. The distribution D =< Xi,X2 > is bracket generating and 
has not a constant rank on R^. If we take the regular Lagrangian 



on E = D and use a result from f83j we obtain a regular Hamiltonian H on 
E* 



2 



2 (l-e2)2 (l-£2) l_^2 

and from (2.3.6) we get 

^''^ = ( ^ \ / Pi 

fl2 J \0 X^ J \p2 

and the corresponding Hamiltonian H on T*M is 

^ ^ 1 / / (pi)2 ^ (P2)2(xi)2 epi ^ 



2 \^\/ (l-e2)2 (1-^2) 1_£2 
From the Hamilton equations (2.3.7) we obtain 



.i_(l + e)^Pi e / (pi)2 a2(xi)2 epf 



(l_e2)2 l_£2y(l_^2)2 l_g2 (1_£2)3 / (piP . a^(xTj2 

Y (1— e-^)-^ 

x^=^-^-'-;^^^^^= (2.3.8) 
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Pr- 



+ 



1 



(l-e2)2 



(1-^ + l-£2 

We make the following change of variables 

' 3.1 ^ vaH!r(t)sinA0(t), 



a2(xi)2 



P2 = a = Ct. 



^ = {I - e^)r{t) COS AO {t), 



and from (2.3.8) we get 



A{1 - e^){^/l^ dO 



dt 



[I - e cos AOf, 



(1 -e2)(7r^dr 



dt 



Hence 



and 



t 



c(l - e cos ^6*) 
A{l-e^)Vl^ 



er sin AO {e cos AO — 1). 
c G i? 
(i0 



^1 -e cos ^61)^ 

But the geodesies are parameterized by arclength, that corresponds to fix 
the level 1/2 of the Hamiltonian and we have 

^2 



H {I - ecos AOf = 



so c = ±1 and therefore 



1 



and finally we obtain 



= ±Vl -e 



1 — e cos AO 
r sin AO 



X 



A{l-e^)Vl^ 



a(l — ecos AO) ' 
sin^ AO 



(2.3.9) 



^1 - ecos ^6*)' 



.dO. 



Remark 2.3.3 l.The above geodesies are the only minimizers of this sub- 
Finslerian geometry. 

2. For e = we obtain the geodesies of the Grushin plane endowed with 
the standard Euclidean metric, i.e. a sub-Riemannian geometry (distribu- 
tional systems with quadratic cost) 

t sin 2at 



, sin at 



(2.3.10) 
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2.3.2 Heisenberg case 



We consider the following distributional system with positive homogeneous 
cost (Heisenberg group) 



X = U^Xi + u'^X2, X 



mm 
«(•) Jo 



J='{u{t))dt, T{u) = \\u\\ + {h, u),b= {e, 0)*, u = (u\ u^)*, < e < 1 



We are looking for the geodesies starting from the origin and parameterized 
by arclength. Here 

d x"^ d 

Xi 



Xo 



dx^ 2 

d x^ d 
+ 



dx"^ 2 dx^ ' 
d 



^3 — [Xl,X2] 



dx'^ 



and {Xi, X2,X^) = TM^, hence the distribution D =< Xi,X2 > of constant 
rank is strong bracket generating [118j . 

Remark 2.3.4 We can try to work directly on the cotangent bundle by com- 
puting the Hamiltonian H{x,p) = , n = cr*{p). Since 



^1 



1 -x72 
1 x^/2 




we obtain 



H 



\ 



1-e 



2\2 



V 



(2.3.11) 



(2.3.12) 



Unfortunately, with this Hamiltonian (2.3.7) is a very complicated system 
of implicit differential equations. 

We will use a different approach. Let us take M = and E = TM 
with the basis {Xi^X2,X3}. E is a Lie algebroid over M with at least one 
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structural function nonzero. The anchor a : E ^ TM is just the identity 
and the matrix of a with respect to the basis of E and TM basis is 




< = 1 I (2.3.13) 

\ 2 2 

Now the above control system can be written 



X = U^Xi +u2X2 + 0X3, X=\ x2 I Gto3 



X 



r r 

min / T{u{t))dt = min / C{u{t))dt, 
"(•) Jo H-) Jo 



where = \J (ui)^ + (n^)^ + eu^ and £ = ■ To solve this problem we 
form the augmented Lagrangian 

£ = /:(n) + An^ 

(A is Lagrange a multiplier) and still work via the maximum principle but 
at the level of the Lie algebroid. If we set 

5 £ ■ ~ 

we obtain 
and therefore 

■H = /iiu*- £= TTT^^ + TT^^^ + An^ - £ - An^ = £, 

because £ is 2-homogeneous. Again, using a result from [83j , the Hamilto- 
nian on E* is given by 

/ \ X 2 



^ = 2 I^V (1372)2 + (13^2)- 13^2 j ' /^3 = A. (2.3.14) 

From the structure equations of Lie algebroids (1.2.8) and the relation 
[Xq,,X/3] = L2/3^7 w6 obtain the non-zero structural functions 

-^12 — 1) -^21 = 
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Now, using the Hamilton equations on Lie algebroids (2.1.3) we get the 
following systems of differential equations 

x^ = |^, x^ = |^, (2.3.15) 



X 



3 



1-2 •12 

XX XX 



and 



2 dm 2 dm 2 



n dn an 
dm dm 

.3 dH . 



where 



A2 = -L^i^s^ = A— (2.3.16) 

dm dm 
m = ^ X = ct. 

dH + e') m + e/x? 



(1-sT (l-.^)V(fe^ + fe? 

/i2 eAtiAt2 



We may use to the following transformations 

^i(t) = (1 - e2)r(t)(acos A^(t) - 6sin^0(t)) 



(2.3.17) 



mit) = Vl -e'^rit){a sin A9{t) + bcosAe{t)) 
such that + 6^ = 1 and we get 



(Mi)^ ^ (/X2)2 



Y (1 - £2)2 (1 _ £2) 

Also, from (2.3.16) we obtain 



^ _ (1 - e'^)VT^A f de 



A J {1- e {a COS AO -bsm. AO)) 

and 

1 



2 ' 



c(l - £(acos^6' - 6sinA6'))' 

But the geodesies are parameterized by arclength, that corresponds to fix 
the level 1/2 of the Hamiltonian and we have 

'H = — (1- e{acosAe -bsmAe)f = 
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so c = ±1 and therefore 

1 



From (2.3.17) we obtain 



^1 - e{acosA9 - bsinA9)' 



(l-e2)(acos A9-bsmAe) 
l-e{acosAe -bsinAe) ' 



But 



so 



Vl - e^(a sin A6i + 6 cos AO) 
1 - e {acosA9 - bsinAO) 

X (c'j = ai. 



A 

Since we are looking for geodesies starting from the origin, we have oi 
and therefore 

X{l—sa) 



x\9) 

and from 



Vl - (a sin A9 + b cos A9) ^\-e^b 
A (1-e (a cos ^6* -6 sin ^61)) ~ A(l - ea) ' 



/ii = — Ax^, 



we obtain 



2. . _ (l-£^) (6 sin ^61-0 COS A(9) (1 - £^)a 
^ ^ ~ A(l-e(acosAe-6sin^0)) ^ A(l - ea) ' 



Finahy, from 



X 



3 



1-2 -12 



2 2 

by straightforward computation we get 

3 , _ (1 -£2)^/1^^^ |- l-cos^6' 



a; (^) = ^ ^ / ^d9. 

2A2(l-ea) ] (\ - e{azo^ A9 -b^ixiA9)f 

Remark 2.3.5 For e = ^ we obtain the sub-Riemannian case (distribu- 
tional systems with quadratic cost) with the solutions 

a sin At — 6(1 — cos X9) 



x{t) 

y{0) 



A 

b sin At + a(l — cos Xt) 
A 



At — sin At n ,n 
40) = ^2 , a^ + b^ = l 



(see also {T^ \3^). 
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